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c3 : Abstract 

We present an analytical calculation of the two-loop QCD corrections to the electromagnetic form 
factor of heavy quarks. The two-loop contributions to the form factor are reduced to linear com- 
binations of master integrals, which are computed through higher orders in the parameter of di- 
mensional regularization e = (4 - D)/2. Our result includes all terms of order e at two loops and 
extends the previous literature. We apply the exponentiation of the heavy-quark form factor to 
derive new improved three-loop expansions in the high-energy limit. We also discuss the implica- 
tions for predictions of massive n-parton amplitudes based on massless results in the limit, where 
the quark mass is small compared to all kinematical invariants. 



1 Introduction 



The electromagnetic form factor of quarks comprises the simplest example of a scattering ampli- 
tude in Quantum Chromodynamics (QCD) and is of basic importance both for theory and phe- 
nomenology. To mention one example, the heavy-quark form factor is directly linked to phe- 
nomenological predictions for the forward-backward asymmetry for inclusive bottom quark pro- 
duction in electron-positron annihilation. In a wider theoretical context, the heavy-quark form 
factor is a quantity of interest in its own. Being gauge invariant, it provides the simplest instance 
for the study of mass effects in QCD hard scattering processes including higher-order quantum 
corrections and, as a key building block, it also relates to many other hard scattering processes 
involving massive quarks. For these reasons the form factor of massive quarks and in particular its 
QCD radiative corrections have received much attention in recent years. 

In an impressive series of papers the two-loop QCD corrections to the heavy-quark form fac- 
tors have been obtained for the vector- and axial-vector coupling as well as the anomaly contribu- 
tions |[Ill2l[3]|. Besides being mandatory for precision phenomenology, these explicit results for 
the massive form factor also exhibit certain universal features of higher order radiative corrections. 
It has long been known, that the leading logarithms of Sudakov type [4] in heavy-quark mass m 
can be resummed [5J, while the exponentiation of the complete coUinear and infrared singulari- 
ties of the massive form factor within dimensional regularization has been addressed only rather 
recently |[6l. 

As an immediate consequence it is now possible to derive important partial information about 
massive n-parton gauge amplitudes from purely massless calculations at any order in perturbation 
theory [6], see also [171 [H. This exploits the factorization properties in the soft and coUinear limit 
of the respective amplitudes, i.e. the one with all partons massless [l9l[T0l| and the amplitude with, 
say, fermions of mass m in the presence of a typical large kinematical invariant in the limit 

10^ 0. The procedure results in a universal and process independent multiplicative relation 
between the two amplitudes. The concept has been put to test in various applications up to two 
loops, e.g. the derivation of the QCD corrections to hadro-production of top-quarks [[m[T2| or the 
Quantum Electrodynamics (QED) correction to Bhabha scattering [8J. Currently, this is an active 
area of research |[T3l[T4ll . 

In the present paper we study the two-loop QCD corrections for the vector form factor of 
heavy quarks in dimensional regularization. We first perform an independent check of Ref. dUl. 
Subsequently, we provide new results by extending the two-loop results in D = 4 - 2e dimensions 
up to terms of order e. To that end the relevant master integrals at two loops have been computed to 
sufficiently high orders in e and the results for those are given as well. We are then in the position 
to discuss several applications. After testing the D-dimensional exponentiation to order e at two 
loops, we can utilize its predictive power to derive a new improved three-loop prediction for the 
form factor to logarithmic accuracy in the heavy-quark mass m. In the same spirit, we extend the 
results for the universal Z factor relating massless and massive amplitudes in the small-mass limit. 
Finally, as a by-product of our calculation, we also investigate the threshold limit of the heavy- 
quark form factor up to two loops and we discuss the Coulomb singularities for color-singlet and 
octet currents. 

The outline of the paper is as follows. In Section [2] we recall the definition of the vector form 
factor of massive quarks and sketch the important parts of our computation. We also provide some 
details of the integral reduction and the relevant master integrals. Our main result, the two-loop 
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QCD contributions for the form factor up to order e, is presented in Sec. [3] along with expansions 
in kinematical limits (threshold and high energy production). Finally, Sec. |4] discusses the above 
mentioned applications, i.e. the exponentiation which leads to a new three-loop prediction for the 
form factor and to improved formulae for relating massless and massive amplitudes through three 
loops. We conclude in Sec. [51 



2 Calculation 

The coupling of a vector current to a heavy-quark pair can be conveniently written in terms of a 
vertex function F^. Given a photon of virtuality (we take space-like = -Q^ < throughout 
this Section) the expression for F^ reads 

^n(h,k2) = ieq<A(^l)|r^^^l(e^m^a,)-^c7■^v^^:^2(e^m^a,)j(A(fc2). (1) 

Here the external quark (anti-quark) of incoming momentum ki (k2) is on-shell with m denoting 
its mass and eq its charge, thus = (and k^ = m^), and it is ^-l-^i -1-^2 = 0. The scalar functions 
and on the right-hand side are the electric and magnetic space-like quark form factors, i.e. 
our main quantities of interest. 




Tree One 



Figure 1 : Feynman diagrams contributing to the quark form factor of massive quarks up to one loop. 

With the help of suitable projections of F^ (see e.g. Ref. [HI) 'Fi and T2 in Eq. ([T]) can be 
computed. They are gauge invariant, but in general, at higher orders in perturbative QCD divergent. 
Both, fi and enjoy a power expansion in the strong coupling as at the scale /i^ (to be employed 
also for all other quantities throughout this article) 

= E^^''-E(''V))v«. (2) 

(=0 ^ ' i=0 

where we have introduced the shorthand notation as = as/{4-n). All necessary Feynman diagrams 
entering the computation of the electric and magnetic form factors fi and fi through NNLO in 
QCD are displayed in Figs. \T\ and [2l At the two loop level, we have summarized the contribu- 
tions with a self-energy correction to the gluon line into the two diagrams QL and GL as follows. 
Diagram QL represents both the heavy and the n, light quarks in the loop correction to the gluon 
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L C QV S 




QL GL GV 

Figure 2: Feynman diagrams contributing to the quark form factor of massive quarks at two loops. The 
diagrams GV, QV and S carry a symmetry factor of 2. Diagram QL denotes both heavy and light quarks in 
the loop and diagram GL represents the gluon and the ghost loop contribution to the gluon self-energy. 



propagator. Likewise diagram GL stands for the gluon and the ghost loop. Throughout this arti- 
cle, = + 1 is the total number of flavors, which is the sum of one heavy {rij^ = 1) and light 
quarks. The notation in Fig.[2]is taken over from the corresponding calculation for the form factor 
of massless quarks, see e.g. Ref. [15J . 



2.1 Master integrals 

The number of Feynman diagrams in Figs. [Hand [2] is relatively small. Thus, their treatment and the 
reduction of the corresponding Feynman integrals proceeds in a standard way. The reduction has 
been achieved with the Laporta-Remiddi algorithm [fT6l[T7l . as implemented in the Maple program 
AIR IfTSll . In this way a list of so-called master integrals, possibly with additional irreducible scalar 
products of loop momenta in the numerator (in short: 'numerators') has been obtained together 
with the corresponding algebraic relations between the master integrals and the form factors T\ 
and Ti- 

We display all master integrals in Figs. [3] and |4] and tabulate them in Tab. [U We found it 
convenient to replace master integrals with numerators by so-called dotted ones. Dotted master 
integrals are free of numerators but have higher powers of propagators. They are independent 
of the (arbitrary) momentum configurations inside the loop graphs. This transformation can be 
achieved with appropriate integration-by-parts (IBP) relations. 

The Feynman integrals are defined as follows: 



^ {ind/2f J (q2 _ . . . (^2 _ ^2yj _ _ _ (^2^ _ ^2^y^ ' 

where L is the number of loops and X some numerator structure. With this definition, for example. 
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Figure 3: The tadpole and the two-point master integrals. The suffix OS denotes on-mass-shell master 
integrals. 
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Figure 4: The three-point master integrals. 



the one-loop tadpole Tlllm becomes: 



Tlllm 



J k^-m^ 



- M) 



1-e 



e'^^ r(-l+6) 



+ 0{e\ 



(4) 



Actually, all master integrals have already been considered in a project concerned with the two- 
loop QED corrections to Bhabha scattering |fT9ll20l . as well as during the previous two-loop QCD 
computation of the heavy-quark form factor in Ref. [T| (see also the References in the last column 
of Tab. [T])- Because we are aiming at the determination of and T2 including the 0{e) terms at 
two loops, the master integrals are needed here to higher powers in e. Therefore we will point out 
specifically the required extensions (see Tab.[I]). 
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Table 1 : The master integrals entering the two-loop vertex diagrams in Fig. [2l The notations follows 
Refs. ITHIlOl. Stars '*' denote one-loop master integrals, the suffix 'OS' on-mass-sheU functions, and the 
suffix 'd' dotted masters. 



For instance, the on-shell master integral SE313mOS (see Fig. [3]) has to be evaluated until terms 
of the order 0{e^). The solution for SE313mOS, which is, together with the other on-mass-shell 
master SESllmOS, actually one of the simplest masters, has been obtained using the package 
0N-SHELL2 [EH, 



SE313mOS 
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4e 
+ 28^3 + 
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596 
12 



36 
5 



+ € 



16 

114181 212196 936 4 

+ 8056 - 9306 + ^ + 62461n2 + 521n^ 2 



32 



32 



128 
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The computation of QCD corrections to form factors is a one-scale problem with variable Q^/nP' 
and the results can be expressed in terms of harmonic polylogarithms (HPLs) up to weight five 
Hmi,...,m,y{^), = 0, + 1. The HPLs were introduced in [i22il . They depend on the conformal 
variable x. 



X = 



(6) 



with 2^ >OandO<;c< 1. 
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One- scale master integrals can be naturally solved by the method of differential equations ||29l 
l30 l which allows the determination of the master integrals analytically to the desired order in e. 
The algebraic relations between masters with numerators and dotted master integrals sometimes 
contain coefficients with additional singularities in e. This is a well-known feature in the algebraic 
reductions based on integrations by part requiring a deeper expansion of the corresponding master 
integrals in e. The complete set of master integrals needed, some of them being quite lengthy, 
contains HPLs up to weight six. We refrain from printing the explicit results in this paper. They 
can be found at the webpage [31 J in the file Masters_V_QCD_21oops_eps .m, 

Let us briefly describe the checks on the master integrals. As was mentioned, we used differ- 
ential equations to compute the analytic .^-dependence of all coefficients of the Laurent expansions 
in 6 in terms of HPLs. Numerical evaluations of the HPLs may be done with the Mathematica 
library HPL [32J. In addition, we evaluated the master integrals at selected kinematical points in the 
Euclidean region by two independent methods. To that end we (i) used the sector decomposition 
package sector_decoinposition [|33l and (ii) derived and evaluated Mellin-Bames representa- 
tions for the masters with the programs ambre [|34l and MB |[35ll . 

In principle, the method of differential equations allows the computation of master integrals to 
any depth in e. In practice, though, every additional power of e poses a new challenge. Neverthe- 
less, we have been able to push the calculation for each of the master integrals even to one power in 
6 beyond what is actually needed (see Tab[T]). For instance, the on-shell master integral SE313inOS 
(see Fig. |3]) has been evaluated until terms of the order 0{e'^). Due to the iterative nature of the 
method of differential equations this implies further constraints and an additional self-consistency 
check on the correctness of the solutions. As a consequence of our strategy we met at intermedi- 
ate steps HPLs up to weight seven. We observe the cancellation of all HPLs with weights bigger 
than five at order 0(e) after summation of all two-loop Feynman diagrams in Fig. [21 An analo- 
gous statement applies to the corresponding constants of weight six. All terms like e.g. ^5, , or 
se = 5(-5,-l;oo) = 0.9874414264032997137... (the latter one is defined e.g. in Ref. (3T\) disap- 
pear in the final result. The master integrals and their Laurent series are now known sufficiently 
deep in € for a computation of the QCD form factor at two loops including all terms of order e^. 
This is clearly beyond the scope of the present paper, but we comment on this issue again in Sec.lH 

2.2 Renormalization 

From the Feynman diagrams of Figs. [U and [2] we obtain the bare results for and within 
dimensional regularization, D = 4-2e. The ultraviolet divergences require renormalization by 
adding the appropriate counter-terms CTi and CT2 and we briefly present all necessary formu- 
lae. The relevant parameters to be renormalized are the strong coupling constant a^, the external 
(heavy-quark) wave function ifr and the heavy-quark mass m. The latter denotes the pole mass, so 
that the renormalization of m imposes the on-shell condition. The bare (unrenormalized) quantities 
are multiplicatively related to the renormalized ones, 

a^,=Za,as, m°=Zinm, if/" = ^fz^ij/ , (7) 

and we work in the MS-scheme, but Z2 and are to be taken in the on-shell scheme. The nec- 
essary renormalization constants are all known [|Tl|36l[37l and through two loops also to sufficient 
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depth in e, 

ZS = (8) 

z;.. . r,L,(gf{-c,(J.^)), 

4" ^ r(i..,(£:r{-c,(^^)), 



7(2) 

^2 



7/y" \ fl /9 11 \ 1 /51 127 11 \ 

,/433 \ / 1705 49 \ 

+ - 24^3 - 78^2 + 96^2ln2j + CpC^(-— + 12^3 + y ^2 - 48f2ln2j 

/I 13 \ /947 \ 

+ CpT^n,[— + 10^2) + CpT^[^ - 30^2) 

r /211 A , 1008 .\ 

+ eCp Cpi — - 2676 - 2946 + 55261n2 - 384Li4(l /2) - 161n^ 2 - 19261n2 2 + —^6 
L \ 16 5 / 

/ 9907 224 398 504 \ 

+ ^ T^^ ^ T^' " 27661n2 + 192Li4(l/2) + 81n4 2 + 9661n2 2 - — 6') 

/851 116 44 \ ^ / 17971 448 340 , \-\] 

where we have always set the factor 5^ = (AnY exp{-eyE) = 1 and in the MS-scheme fio = 
ll/3C^-2/3n^. 

In addition to Eqs. (f8l)- (fTT]) the two-loop counter-terms CTj^\ i = 1,2 contain the subtraction 
of those one-loop sub-divergences from the two-loop diagrams of Fig. [2] which are related to the 
on-shell renormalization of the heavy mass m in the loop. At one loop we denote the relevant 
results by f^'^^'^^'' and f^'^^'^^^ They can be derived from diagram One in Fig.[I]upon insertion of 
a one-loop mass counter-term for the heavy quark in the loop (see 111), f^'^^'^^'' and need 
to be calculated to order O(e^) and we find in space-like kinematics, = -Q^ < 0, the following 
results 



y^j"""'^^^ = -Cf](l-2(;c3)-F3;c2-3xi)(2//o) + 2//o7?W-2(l-F2;c2-2xi)[ (12) 
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+ Cp\4 + (-26 - 4//_i,o + 2Ho + lHofl)p{x) + (-26 - 4//_i,o + 2//o,o)(l - 2{x^) + 3x2 - 3xi ) 
+ 27/o(l-4(x3) + 6;c2-4;ci)l 



+ eCpl^4Hoil - 6(X3) + 9x2 - 5xi) - 8(1 + 3x2 - 3xi) + (-46 - 26 + 4//„i6 + 8//-i,-i,o 
- 4//-i,o - 4//_i,o,o + 4Hq - //06 - 4//o,-i,o + 2//o,o + 2Ho,iio)pix) + (-46 + 4//_i6 
+ 8//-i,-i,o - 4//_i,o,o - Ho(2 - 4//o,-i,o + 2//o,o,o)(l - 2(x3) + 3x2 - 3xi) - 36(1 - 8/3(x3) 

+ 14/3x2 - 10/3x1 ) - (4//_i,o - 2//o,o)(l -4(x3) + 6x2 - 4xi)l 



+ 6"Cf jl6(l + X2 - xi) + (-46 - 46 - 14/56^ + 8//-16 + 4//-16 - 8//-i,-i6 - 16//-i,„i,-i,o 
+ 8//- 1-1,0 + 8//_ 1-1,0,0 - 8//_i,o + 2//-i,o6 + 8//-i,o,-i,o -4//_i, 0,0 -4//_ 1,0,0,0 + 8//0 
- I4/3//06 --^^06 + 4//o,-i^2 + 8//o,-i,-i,o-4//o,-i,o-4//o,_i,o,o + 4//o,o -^0,06 -4//o,o,-i,o 
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^mct,(i) ^ -CpU{x2-xi) + 4HQ{x^-3l2x2 + ll2xi)-HQq{x)\ (13) 



+ 2//0AO + 2//o,o,o,o);?(x) -2^2(1 - 12(X3) + 18x2 - lOxi) - 10/3^3(1 -24/5(x3) + 34/5x2 
- 22/5xi) + (-14/5^2' + - 8//-i,-i6 - 16//-i,„i,^i,o + 8//-i,.i,o,o + 2//.i,o^2 

+ 8//-i,() -1,0 -4//_ 1,0,0,0- 14/3//o^3 +4//o-i^2 + 8//o,-i,-i,o-4//o,-i,o,o -//o,o^2 
-4//o,o,_i,o + 2//o,o,o,o)(l - 2(X3) + 3x2 - 3xi) + (4//_i^2 + 8//_i,_i,o -4//_i,o,o + 8//0 -^06 

- 4//o,-i,o + 2//o,o,o)(l - 4(X3) + 6x2 - 4xi) - (8//^i,o - 4//o,o)(l - 6(x3) + 9x2 - 5xi )| + 0{e^) , 
1 
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+ (8(X2 - xi) + (I6//0 + 4//o,o - 4^2 - 8//-i,o)(x3 - 3/2x2 + 1/2x0 + iz2 + 2//_i,o - //o,o)^(x)J 

+ eCf {8(X2 - xi) + (8//_i^2 - 8^3 + 16//-i,-i,o - 32//_i,o - 8//-i,o,o + 24Ho - 2Ho^2 - 8//o,-i,o 
+ 16//o,o + 4//o,o,o)(x3 - 3/2x2 + 1/2x0 + (2^3 - 2//-i^2 - 4//-i,-i,o + 2//_i,o,o + l/2//o^2 
+ 2//o,-i,o - Ho,o,o)qix) - 16^2(^3 - 13/8x2 + 5/8xol 



+ e^Cf |l6(x2 - xO - 24^2(-^3 - 5/3x2 + 2/3x0 - 32^3(x3 - 35/24x2 + 1 l/24xi ) 

+ (-28/5^2' + 16//-i^3 + 32//_i^2 - 16//-i,-i^2 - 32//_i,_i,_i,o + 64//_i,_i,o + 16//_i,_i,o,o 
- 48//-i,o + 4//_i,o^2 + 16//-i,o,-i,o - 32//_i,o,o - 8//-i,o,o,o + 32//o - 28/3//o^3 - 8//o^2 
+ 8//o,-i^2 + I6//0,- 1,-1,0 - 32//o,-i,o - 8//0,- 1,0,0 + 24//o,o - 2//o,o^2 - 8//o,o,-i,o + 16//o,o,o 
+ 4//o,o,o,o)(^3 - 3/2x2 + 1/2x0 + (V/5^2^ - 4//-i^3 + 4//-i,-i6 + 8//_i,_i,_i,o -4//_i,_ 1,0,0 
-^-1,0^2 -4//-i,o,-i,o + 2//-i,o,o,o-4//o + 7/3//o^3 -2//o,-i^2 -4//o,- 1,-1,0 + 2//o,- 1,0,0 

+ l/2//o,o^2 + 2//o,o,-i,o - //o,o,o,o)'?(^)} + O(e^) ■ 

^mct,(i) ^mct,(i) computed before to order 0(6) (see flU), the 0(6^) terms are newQ 

All results are given in terms of HPLs with argument < ,x; < 1 . Moreover, it is convenient to 
abbreviate the ubiquitous polynomial 

\ + 1 

= = 2(^qqW+;?qq(-^)), (14) 

which is related to the quark-quark splitting function at one-loop, Pqqix) = 2/(1 - x) - 1 - We 
use also the additional definitions 

(1 +X)" 

q(x) = (16) 

I -X l+x 

With these ingredients we arrive at the following results for counter-terms at one- and two- 
loops, CT^^^ and CTf'\ These counter-terms need to be added to the sum of the unsubtracted 
(bare) diagrams in order to arrive at the (ultraviolet) renormalized heavy-quark form factors. 

CTf^ = zfTf\ (17) 



'There are some misprints in Eqs. (37) and (38) of Ref. ( 1]. However the source files of the authors agree with the 
corresponding parts of our Eqs. (fT2] i and ( fTST l. 
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where / = 1,2 and, trivially, T'^^ = 1 and = at the Born level. The bare results at one loop 
are needed to order 0{e^) as well. We do not display them here, but the interested reader can 
easily derive them with the help of Eqs. (fTTl) and (fT9l ) from the renormalized result T'^^^ presented 
in the next Section. To order 0{e) the explicit expressions for f^'''^^ and are also given in [1 J. 

3 Results 

We are now in a position to present our main results. In space-like kinematics, = < 0, the 
renormalized form factors 'J^i (electric) and (magnetic) through NNLO in QCD and up to 0{e) 
are expressed in terms of HPLs of argument x, < x < \, see Eq. We also factorize again 
the polynomials p{x), q{x) and employ the abbreviations jc„, see Eqs. (fT4l) - (fT6l) . Adding the bare 
diagrams in Figs. [Hand [Altogether with the counter-terms of Eqs. (flTl) and (fTSi) . 

<jr{l) ^ ^^l^'W+Crf*, (19) 

and setting the scale to /i^ = rrp-, we have: 

rl^^ = (20) 

-^Cp\-2-2p{x)H^\ 

- Cp\a + (1 - 2xi)Ho + 2p{x){Hofl - 2//_i,o - ^2 + 2//o)J 

- eCp^% + p{x)mo + 8//-i,_i,o -4^3 +4//_i^2 -4//-i,o.o -^0^2 -4//0.-1.0 + 2//o,o,o + 4//o,o 
- 8//_i,o - 4^2) + (1 - 2xi)(//o,o - 2//-i,o) + 2x1^2) 



^Cp^\6 + 2^2 - 2/3^3 + P{x){-H2 - 86 - 7^ + 8//_i6 - 8//_i,_i6 - 16//_i,_i,_i,o + 8//-i,-i,o,o 
+ 8//-i,o,-i,o + 2//-i,o6 -4//_ij),o,o+ I6//0 + 8//-16 + 16//_i _i_o - 8//_i,o,o -2//o6 - 8//()-i,o 
+ 4//o,o,o- I4/3//06 +4//o,-i6 + 8//o,-i,-i,o-4//o,- 1,0,0- 16//-i,o - ^0,06 -4//o,o,- 1,0 + 2//o,o,o,o 
+ 8//o,o) + (1 - 2xi)(2//_i6 + 4//_i,_i,o - 2//_i,o,o - I/2//06 - 2^ - 2//o,-i,o + «o,o,o)} + O(e^) , 



^(2) 
^ 1 



^C/{ 2 + 4p(x)//o + 4p(x)2//o,o 
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+ ^CfC^{ll/3 + ll/3Mx)//o 
- ^Cf?i,r/,{4/3 + 4/3K^)//o 
+ -CpA% + 4p(x)(4//o - 6 - 2//-i,o + //o,o) + 2(1 - 2xi)(//o - 2//o,o) - 8x2//o,o 



+ Ap{xf{AH^,fl + 3//o,o,o - 6^0 - 4//^i,o,o - 2//o,-i,o) 
+ C^{-49/9 + 4^ + 4//o,o + p{x)(4(2Ho - 67/9//o - 2^ + 4//_i,o - 4//o,o - 4//i,o + 4//o,o,o) 



+ Ipixf i2Ho,-i,o - ^3 - (iHo - 2//o,i,o - 2//o,o,o) j 
^Cf«,r^ {20/9 + 20/9p(x)Ho^ 

CpTjfi^{\ - 24x5 + 60x4 - 44x3 + 6x2)(4//o^2 + 4//o,o,o) - 8/3^2(1 - 196x4 + 392x3 - 229x2 
+ 33xi) + 76/9//o,o(l +248/19x4-496/19x3 + 326/19x2 -78/19xi) + 46/3//o(l +712/69x3 

- 356/23x2 + 218/69x1 ) + 1532/27 + 784/9x2 - 784/9xi + (944/27i7o + 8//06 + 20/3Ho,o,o)p(x)^ 

C^C^|-1595/27 + 12x2 - 12xi +(1 -4688/37x5 + 11720/37x4-29924/111x3+3242/37x2 

- 1 120/ lllxi)(lll/40^2^) + (l -2256/29x5 +5640/29x4 -4684/29x3 + 1386/29x2 

- 144/29xi)(29/2//o -1,0,0) + (1 -648/11x5 + 1620/11x4 - 1324/11x3 + 366/11x2 -36/llxi)(-227^o^3) 
+ (1 - 56x5 + 140x4 - 344/3x3 + 32x2 - 10/3xi)(24//i,o^2 + 24//i,o,o,o) + (1 - 368/7x5 + 920/7x4 

- 2252/21x3 + 206/7x2 - 64/21xi)(21/4//o,o6) + (1 - 240/7x5 + 600/7x4 - 484/7x3 + 18x2 

- 16/7x1 )(21/2//o,-i ^2) + (1-516/17x5 + 1488/17x4 - 1558/17x3 + 42x2 - 156/17xi)(34i7o6) 
+ (1 - 24x5 + 60x4 - 44x3 + 6x2)(-8//o,o,i,o) + (1 - 172/9x5 + 544/9x4 - 1846/27x3 + 872/27x2 

- 170/27xi)(5477o,o,o) + (l +48/13x5 - 120/13x4 + 196/13x3- 174/13x2 + 24/13xi)(1377o,o,-i,o) 

+ (1 + 8x5 - 20x4 + 24x3 - 16x2 + 2xi)(12i7o,i,o,o) + (1 +48x5 - 120x4 + 68x3 + 18x2 - 16xi)(l/477o,o,o,o) 

+ (1 - 288/89x4 + 468/89x3 - 18/89x2 - 190/89xi )(-89/377o,o) - (1 - 3x4 + 6x3 - 19/4x2 + 7/4xi)(3277i,o,o) 

+ (1 - 12/1 1x4 + 24/1 1x3 - 14/1 1x2 + 2/1 lxi)(-44i7o,-i,o) + (1 + 12x4 - 24x3 + 14x2 - 2xl)(16i/o,i,o) 

+ (1 + 15x4 - 30x3 +45/2x2 - 15/2xi)(-24//_i^2) - (1 + 32x4-484/7x3 + 362/7x2-926/63x0(21^2) 

+ (1 + 47x4 - 94x3 + 1 13/2x2 - 19/2xi )(-24//_i,o,o) + (1 + 54x4 - IO8X3 + 187/3x2 - 25/3x0(24^3) 

+ (1 - 24x3 + 36x2 - 14xi)(-4//i,o) + (1 - 648/125x3 + 972/125x2 - 574/125xl)(-125/18//o) 

+ (1 + 9/2x3 - 27/4x2 + l/4xi)(16/3//-i,o) + (1 + 6x2 - 6xi)(24^2ln2) + (62/3^3 + 158/9^2 + 383/40^2^ 

- 32/3//_i^2 - 208/3//_i,_i,o + 172/9//_i,o + 8//-i,o6 + 140/3//-i,o,o + 8//-i,o,o,o + 24//_i,i,o 

- 1460/27//0 + 6//o^3 +44/3//o^2 -27/2//o,-i^2 + 140/3//o,-i,o - 19/2//o,-i,o,o -50/9//o,o + 21/4//o,()6 

- 27//o,o,-i,o - 16/3//o,o,o + 97/4//o,o,o,o + 16//oAi,o - 16//o,i,o + 4//o,i,o,o + 4//i^2 + 24//i,_i,o + 8//i,o 

- 16i7i,o,o - 8i7i,i,o)p(x) + (-37/5^2^ + 877_i6 + 8//-i,o^2 - 16i/-i,o,-i,o + 16i/-i,oAO + 16i/-i,o,i,o - 2i/o6 
+ 4i/o,-i^2 - 40/^0,-1,-1,0 + 32/^0,-1,0,0 + 24//o,-i,i,o + 4//o,o^2 + 48//o,o,-i,o - 24//o,o,o,o - 32//o,o,i,o + 4^o,i^2 

+ 24//o,i,-i,o - 16//o,i,o,o - 8//o,i,i,o - 8//i^3 - m,o(2 + 16//i,o,-i,o - 16//i,o,o,o - l6Hi^oxoMxf] 

C/|46 + (l - 1008/11x5 +2520/11x4 - 188x3 +582/ 11x2 -48/ 11x0(1 l//o,-i,o,o) + (l - 13680/229x5 

+ 34200/229x4 - 27468/229x3 + 7002/229x2 - 512/229x0(229/20^2^) + (1 - 240/7x5 + 600/7x4 
-484/7x3 + 18x2 - 16/7xi)(-21//o,-i6) + (1 - 768/23x5 + 1920/23x4 - 1492/23x3 + 318/23x2 

-24/23xi)(-46//o,o,-i,o)-(l -360/11x5 +900/11x4 -700/11x3 + 150/11x2 - 12/llxi)(44//o,i,o,o) 
+ (1 -24x5 +60x4-44x3 + 6x2)(16//o,o,i,o) + (l + 120/37x5 -432/37x4+512/37x3 - 170/37x2 

- 88/37x1 )(-74//o,o,o) + (1 + 120/3 Ixj - 3 12/3 IX4 + 272/31x3 - 28/3 1x2 - 1 10/3 Ixi )(-62//o^2) 

+ (1+48x5 - 120x4 + 68x3 + 18x2- 16xi)(-l/2//o,o,o,o) + (l +72x5 - I8OX4 + 156x3 -54x2 + 4xi)(-8//i,o^2 

- 8//i,o,o,o) + (1 + 168x5 - 420x4 + 364x3 - 126x2 + 12xi)(4//o^3) + (1 + 528x5 - 1320x4 + 1 124x3 - 366x2 

+ 32xi)(-l/2/?o,o6) + (1 - 168/5x4 + 336/5x3 - 194/5x2 + 26/5xi)(-20^3) + (1 - 63/2x4 + 63x3 - 127/4x2 
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+ 1/4x1 )(32//_i,o,o) + (1 - 192/13x4 + 21 12/65x3 - 240/13x2 + 14/65xi)(65^2) + (1 + 768/61x4 - 2016/61x3 
+ 1800/61x2 -518/61xi)(61//o,o) + (l + 15x4-30x3+45/2x2 - 15/2xi)(48//_i^2) 

+ (1 + 45/2x4 - 45x3 + 197/8x2 - 17/8xi )(64//i,o,o) + (1 + 192/7x4 - 384/7x3 + 221/7x2 - 29/7xi)(56//o -i,o) 
+ (1 +48x4 - 96x3 + 54x2 - 6xi)(-8//o,i,o) + (1 + 384/17x3 - 576/17x2 + 158/17xi)(34//_i,o) 
+ (1 + 6x2 - 6xi)(-48^2ln2) + (1 - 2xi)(-9/2//o) + - 3/2x2 + l/2xi)(-1927^i,o) + (-16^3 + 122/5^2^ 
+ 6477_i,_i,o,o + 16//-i.o^2 + 32i7_i,o,-i,o - 64//_i,o,o -4877_i,o,o,o - 327^0^2 + 8//o,-i6 + I6//0.-1.-1.0 

- 4077o,-i,o,o + 64//o,o + 2477o,o^2 - 40//o,o,-i,o + 16i/o,o,o + 6877o,o,o,o + 16//o,o,i,o - 32i7o,i,o + 16i/i^3 

+ 16771,0^2 - 32771,0,-1,0 + 32771,0,0,0 + 32Hifixo)p(xf + (-8^ + 8^ - 731/20^2^ + 8i7-i^2 + 16i7-i,-i,o 

+ 14477_i,o - 8i/_i,o,o + 38//0 - 28/?o6 - 2277o6 + 2777o,-i^2 - 8//o,-i,o + 27i/o,-i ,0,0 - 104i7o,o - 65/27?o,o6 

+ 2i?o,o,-i,o- 18/^0,0,0 -81/2/^0,0,0,0 - 16//o,o,i,o-28//o,i,o,o -64//i,o-8i7i,o^2 -8i?i,o,o,o)/'W| 

n;C/,r^{424/27 + (1 - 2xi)(-8/3//_i,o + 2/9//0 + 4/3//o,o) + 4^2(1 + 2/3xi) 

+ (-16/3^3 - 88/9^2 + 16/3//_i^2 + 32/3//_i,_i,o - 176/9//_i,o - 16/3//_i,o,o + 424/27//o 

+ 8/3i/o^2 - 16/3i/o,-i,o + 88/9i7o,o + S/3Ho,o,o)p(x)^ 

6Cf r^{(l - 24x5 + 60x4 - 44x3 + 6x2)(-2/5^2^ + 8//-i,o^2 + 8//_i,o,o,o - SHo^s - 24Ho,-i^2 

- 8//o,-i,o,o + 4//o,o6 - 8//o,o,-i,o + 12//o,o,o,o + 16//o,i,o,o) + (1 - 198/25x5 + 743/25x4 - 1009/25x3 
+ 1201/50x2 -297/50xi)(200/9//o^2) + (l -792/173x5 + 3964/173x4 -6020/173x3 +3706/173x2 
-900/173xi)(346/9i/o,o,o) + (1 -232/7x4 +464/7x3 -262/7x2 + 30/7xi)(224/3^2ln2) 

+ (1-3252/125x4 + 1372/25x3 -4287/125x2 + 122/25xi)(-500/9^2) + (l - 1376/69x4+2752/69x3 

- 1508/69x2 + 44/23x1 )(- 184/3^3) + (1 + 19/3x4-59/7x3 + 163/84x2 -2/3x0(1 12/377o,o) 

+ (1 +248/19x4-496/19x3 +326/19x2 - 78/19xi)(-152/377_i^2 - 152/9i/_i,o,o- 152/9i7o,-i,o 
+ 304/97/1,0,0) + (1 + 11352/1345x3 - 17028/1345x2 + 2986/1345x0(1345 /27//o) + (1 +712/69x3 
- 356/23x2 + 218/69xi)(-92/3//-i,o) + 4138/27 + 1432/9x2 - 1432/9xi + (-944/27^2 - 2^2^ 

- 1888/27//_i,o + 32/3//_i,o^2 + 40/3//_i,o,o,o + 3064/27//o - 128/9//o^3 - 50/9//o^2 - 40//o,-i^2 

- 40/3//o,-i,o,o + 944/27//o,o + 8//o,o6 - 40/3//o,o,-i,o - 74/9//o,o,o + 20//o,o,o,o + 80/3//o,i,o,o);?(x)} 

6Cf C^|(-28745/162 - 1037x6 + 31 1 1x5 - 31 1 1x4 + 1037x3 + 78x2 - 78x0 + (1 - 1200x5 + 3000x4 

- 2572x3 + 858x2 - 88xi)(2//o,i,o,-i,o) - (1 - 976/3x5 + 2440/3x4 - 2084/3x3 + 686/3x2 - 24xi)(9//o,i,o,o,o) 

+ (1 -4688/37x5 + 11720/37x4-29924/111x3+3242/37x2 - 1120/1 11x0(1 ll/20//_i^2^) + (1 - 3664/29x5 
+ 10360/29x4 - 10172/29x3 + 3938/29x2 - 488/29xi)(87/2/ro^3) + (1 - 120x5 + 300x4 - 248x3 + 72x2 - 6x0 
(SHoxoxo) + (1 - 2832/29x5 + 7080/29x4 - 124220/609x3 + 12550/203x2 - 3856/609xi)(-609/40//o^2^) 
-(1 -4240/49x5 + 10936/49x4 - 1412/7x3 +510/7x2 -452/49xi)(98//i,o6) + (1 -9968/121x5 
+ 24920/121x4 - 62452/363x3 + 6306/121x2 - 2048/363x0(1089/8^5) + (1 - 18256/229x5 + 45640/229x4 

- 114628/687x3 + 11674/229x2 -3920/687xi)(-687/20//i ^2^) + (1 -2256/29x5 +5640/29x4 -4684/29x3 
+ 1386/29x2 - 144/29xi)(29//_i,o,-i,o,o)-(l - 1296/17x5+3240/17x4 - 13548/85x3 +4122/85x2 

- 464/85x1 )(85/2//i,o,o^2) + (l - 1964/29x5 + 6368/29x4-7270/29x3+3312/29x2 -432/29xi)(-58//o,o^2) 
+ (1 - 5040/79x5 + 12600/79x4 - 10372/79x3 + 2958/79x2 - 304/79xi)(-79/?o,-i,-i,o,o) 
+ (1 - 312/5x5 + 156x4-652/5x3 + 198/5x2 -24/5xO(-20/?i,o6 -40/?i,o,i,o,o) 
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+ (1 -312/5x5 + 156x4-636/5x3 + 174/5x2 - 16/5xi)(40//o,-i,o,i,o) 

+ (1 -7728/125x5 + 3864/25x4 - 15852/125x3 +4458/125x2 -448/125xi)(125//o,_i,o,o,o) 

+ (1 - 1008/17x5 + 2520/17x4 - 2068/17x3 + 582/17x2 - 60/17xi)(-34//o,i,o6) + (1 - 648/1 1x5 

+ 1620/11x4 - 1324/11x3 + 366/11x2 -36/llxi)(-44//_i,o^3) + (l -4080/71x5 + 10200/71x4-8132/71x3 

+ 1998/71x2 - 128/71xi)(71/4//o,o,o6) + (l -5520/97x5 + 13800/97x4 -11260/97x3+3090/97x2 

- 304/97x1 )(-97//o.-i,o,-i,o) + (1 - 56x5 + 140x4 - 344/3x3 + 32x2 - 10/3xi)(4877_i,i,o6 + 48i7_ 1,1,0,0,0 
+ 48i7i,_i,o^2 + 4877i,_i,o,o,o + 487?i,o,o,i,o -48771,1,0^2 -4877i,i,o,o,o) + (1 -4176/77x5 + 10440/77x4 

- 8556/77x3 + 342/1 1x2 - 256/77xi)(-154//i,o,o,-i,o) + (1 - 368/7x5 + 920/7x4 - 2252/21x3 + 206/7x2 
-64/21xi)(21/2i7_i,o,o6) + (1 -2256/43x5 +5640/43x4 - IO8X3 + 1326/43x2 - 152/43xi)(-43/2i7o,o6) 
+ (1 - 4080/89x5 + 10200/89x4 - 8252/89x3 +2178/89x2 - 224/89xi)(89/2i7o,-i^3) + (1 - 3504/79x5 

+ 8760/79x4 - 6980/79x3 + 1710/79x2 - 144/79xi)(79/2//i,o,o,o,o) - (1 -4240/99x5 + 10480/99x4 

- 8972/99x3 + 2978/99x2 - 316/99xi)(198//i,o,o,o) - (1 -40x5 + 100x4 - 230/3x3 + 15x2 - 1/3x0 
(24//o,OAi,o) + (1 - 648/17x5 + 1620/17x4 - 76x3 + 318/17x2 - 32/17xi)(68//o,-i,o^2) + (1 - 1296/37x5 

+ 3240/37x4 - 2556/37x3 + 594/37x2 - 56/37xi)(74//o,i,i,o,o) + (1 - 240/7x5 + 600/7x4 - 484/7x3 + 18x2 

- 16/7xi)(21//_i,o,_i^2) + (l - 144/5x5 +72x4- 1972/35x3 +438/35x2 -48/35xi)(-105//o,-i,-i^2) 
+ (1 -624/23x5 + 1560/23x4 - 1188/23x3 +222/23x2 - 16/23xi)(-69//i,o,_i^2) + (1 - 336/13x5 

+ 840/13x4-588/13x3 +42/13x2 + 16/13xi)(-1377i,o,-i,o,o) + (1 -24x5 +60x4-228/5x3 +42/5x2 

- 4/5xi)(-20i7o,o,i ,0,0) + (1 - 24x5 + 60x4 - 44x3 + 6x2)(-16//_i,o,o,i,o + 48//o,o,-i,i,o + 8//o,o,i^2 

+ 48770,0,1,-1,0- 16^^0,0,1,1,0) + (1 - 1168/53x5 +4168/53x4 -5132/53x3 +2530/53x2 -436/53x0 
(106//o,o,i,o) + (l -240/11x5 + 600/11x4-532/11x3 + 18x2 -48/1 1x0(1 l/477o,o,o,o,o) 
+ (1 - 1460/67x5 + 4346/67x4 - 4782/67x3 + 33x2 - 405/67xi)(6777o^2) + (1 - 3792/181x5 + 9480/181x4 
-6924/181x3 +906/181x2 -32/181xi)(181/2i7o,o,-i,o,o) + (1 - 144/7x5 + 360/7x4-260/7x3+30/7x2) 
(42//o,-i^2ln2 - 14//o,-i,i,o,o + 42//o,i^2ln2 -42//o,i,-i6 - UHox-mo) + (1 - 3256/175x5 + 9736/175x4 

- 10152/175x3 +4222/175x2 -648/ 175xO(-350//o,_i,o,o) + (l - 1296/77x5 + 3240/77x4-2028/77x3 

- 18/7x2 + 128/77xO(77/2//o,o,o,^i,o) + (l -2944/223x5 +9016/223x4 -9852/223x3 +4622/223x2 

- 1216/223xi)(223/2//o,o,o,o) + (1 - 688/55x5 + 1048/55x4 + 252/55x3 - 826/55x2 + 8/5xO(55//o,i,o,o) 
+ (1 - 104/9x5 +7328/171x4-9440/171x3 + 1654/57x2 -52/9xO(-342//o,o,^i,o) + (l -492/53x5 

+ 1728/53x4-2106/53x3 + 1106/53x2 -332/53xi)(-212//o,-i^2) + (l -464/51x5 + 1160/51x4 

- 1948/153x3-62/17x2 + 112/153xi)(-153/2//o,i^3) + (1 -624/109x5 + 1560/109x4-876/109x3 
-246/109x2 - 32/109x0(109/4^26) -(1 -48/11x5 + 120/11x4 - 116/55x3 -426/55x2 + 72/55x0 
(110//o,o,-i.-i,o) + (l +48/13x5 - 120/13x4 + 196/13x3 - 174/13x2 + 24/13xi)(2677_i,o,o,-i,o) + (l + 8x5 

- 20x4 + 24x3 - 16x2 + 2xi)(24i7_i,o,i,o,o) - (1 +516/23x5 + 306/23x4 - 2030/23x3 + 1479/23x2 - 243/23x0 
(92//_i,o,o,o) + (l + 1392/41x5-3480/41x4 + 3332/41x3 - 1518/41x2 + 192/41xi)(-41/277o,o,-i6) + (l 

+ 2920/61x5 - 9316/61x4 + 10092/61x3 - 3686/61x2 + 92/61xi)(-61/27?o,o,o) + (1 + 48x5 - 120x4 + 68x3 
+ 18x2- 16x0(1 /2//_i,o,o,o,o) + (l + 1117/12x5-2851/12x4 + 25633/120x3 - 1491/20x2 + 659/120x0 
(48^2) + (1 +516/5x5 - 978/5x4 + 538/5x3 - 39/5x2 - 9/5xi)(-20//-i,o6) + (1 - 30x4 + 6OX3 - 38x2 + 8x0 
(80//i,o,_i,o) + (1 - 30x4 + 60x3 - 75/2x2 + 15/2xO(-32//i,o,i,o) + (1 - 736/63x4 + 1268/63x3 - 478/63x2 

- 10/7x1 )(126//o,-i,o) + (1 - 34/3x4 + 44/3x3 - 4xO(-48//o,i,o) + (1 - 2592/467x4 + 5832/467x3 
-4968/467x2 + 1668/467x0(467/3^) + (1 +60/37x4 - 120/37x3+70/37x2 - 10/37xi)(29677o,-i,-i,o) 
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-(1 + 2178/571x4 -7704/571x3 + 18585/1142x2 -4310/571xi)(1142/9//o,o) + (1 + 348/79x4 -696/79x3 
+ 336/79x2 + 12/79xi)(158//i^3) + (l + 196/19x4-488/19x3+402/19x2 - 118/19xi)(-76//i,o,o) 
+ (1 + 12x4 - 24x3 + 32/3x2 + 4/3xi)(-72//^i^2ln2 + 24//_i,i,o,o - 72//i^2ln2 + 72//i _i^2 + 24//i,„ 1,0,0) 
+ (1 + 12x4-24x3 + 14x2 -2xi)(-96//o,-i, 1,0 - 16//o,i^2 -96//o,i -1,0 + 32//o,i,i,o) + (1 + 1542/125x4 

- 3444/125x3 + 5883/250x2 - 3061/375xi)(-250/3^2) + (1 + 14x4-28x3 + 167/9x2 -41/9xi)(216//_i,_i^2) 
+ (1 + 108/5x4 - 216/5x3 + 126/5x2 - 18/5xi)(-120//i, 1,0,0) - (1 + 392/17x4 - 580/17x3 + 86/17x2 

+ 138/ 17x1 )(3477_i,o,o) + (1 + 1020/37x4-2040/37x3 + 1320/37x2 - 300/37xi)(-7477_i^3) 

+ (1+630/17x4 - 1260/17x3+755/17x2 - 125/17xi)(13677_i,_i,o,o) + (l + 115/3x4-230/3x3+835/18x2 

- 145/18xi)(14477_i,o,_i,o) + (1 +312/5x4 - 624/5x3 + 374/5x2 - 62/5xi)(-40//_i,o,i,o) + (1 + 72x4 - 144x3 
+ 132x2 - 60xi)(12^2ln2) + (1 + 296/3x4 - 172x3 + 214/3x2 + 2xi)(-18i7_i^2) + (1 - 22680/473x3 

+ 34020/473x2 - 12286/473xi)(-473/108//o) + (1 - 24x3 + 36x2 - 14xi)(24//_ 1,1,0 + 4//i^2 + 24//i,_i,o 

- 8//i,i,o) + (1 - 164/9x3 + 82/3x2 - 100/9xi )(-36//i,o) + (1 - 7488/763x3 + 1 1232/763x2 - 5270/763x0 
(763/9//_i,o) + (1 - 5x3 + 15/2x2 - 9/2xi)(-48//_i,_i,o) + (1 + 6x2 - 6xi)(-41n'* 2 - 96Li4(l/2) -48^2ln^ 2) 
+ (-137^5 + 22^2^3 + 148/5//_i6^ - 32//_i,_i^3 - 32//_i,_i,o^2 + 64//_i,_i, 0,-1,0 - 64//_i,_ 1,0,0,0 
-64//_i,_i,o,i,o + 8//- 1,0^3 - 16//- 1,0,- 1^2 + 160//- 1,0,- 1,-1,0 - 128//- 1,0,- 1,0,0 -96//- 1,0,- 1,1,0 - 16//- 1,0,0^2 

- 192//_i,o,o,-i,o + 96//- 1,0,0,0,0 + 128//- 1,0,0, 1,0 - 16//-i,o,i^2 - 96//- 1,0, 1,-1,0 + 64//_ 1,0, 1,0,0 + 32//_i,o,i,i,o 

+ 32//_i,i^3 + 32//_i,i,o6 - 64//-i,i,o,-i,o + 64//_ 1,1, 0,0,0 + 64//_ 1,1, 0,1,0 + 24/5//o^2^ - 36//o,-i6 + 8//o,-i,-i6 
+ 304//o,-i,-i,_i,o - 224//o,-i,-i,o,o - 144//o,-i,-i,i,o - 4//o,-i,o^2 - 288//o,-i,o,-i,o + 160//o,-i,o,o,o 
+ 160//o,-i,o,i,o - 24//o,_i,i^2 - 144//o,_i,i,_i,o + 96//o,_i,i,o,o +48//o,_i,i,i,o + 76//o,o^3 - 72//o,o,-i^2 

- 384//o,o,-i,-i,o + 232//o_o__i_o,o + 192//o,o,-i,i,o + 84//o_o,o^2 + 256//o,o,o,_i,o - 48//o,o,o,o,o - 136//o,o,o,i,o 
+ 32//o,o,i^2 + 192//o,o,i,_i,o - 104//o,o,i,o,o - 64//o,o,i,i,o + 12//o,i^3 - 24//o,i,_i^2 - 144//o,i,_i,_i,o 

+ 96//o, 1,-1, 0,0 + 48//o,i -1,1,0 - 20//o,i,o^2 + 160//o,i,o,_i,o - 104//o4,o,o,o - 96//o,i,o,i,o + 8//o,i,i^2 
+ 48//o,i,i,-i,o - 32//o,i,i,o,o - 16//o,i,i,i,o - 148/5//i^2^ + 32//i,_i^3 + 32//i,_i,o^2 - 64//i,_i, 0,-1,0 
+ 64//i,-i,o,o,o + 64//i,_i,o,i,o- 8//i,o^3 + 16//i,o,-i^2 - 160//i,o,-i,-i,o + 128//i,o,-i,o,o + 96//i,o,-i,i,o 
+ 16//i,o,o^2 + 1 92//1, 0,0,- 1, -96//i, 0,0,0,0- 128//i,o,o,i,o + 16//i,o,i^2 + 96//i,o,i,-i,o -64//i,o,i,o,o 

- 32//i,o,i,i,o - 32//i,i^3 - 32//i,i,o^2 + 64//i,i, 0,-1,0 - 64//i,i,o,o,o - 64//i,i,o,i,o)K^)^ 

+ (929/8^5 +356/9^3+2524/27^2 + 13/4^2^3+77^2^ - 156//-i^3 -536/9//-i^2 + 383/20//_i^2^ 
+ 96//_i,-i^2 +480//-i,-i,-i,o -208/9//-i,-i,o + 16//-i,-i,o6 - 312//_i,_i,o,o + 16//-i,-i,o,o,o 

- 144//_i,_i,i,o + 5912/27//-i,o + 12//_i,o6 + 6//_i,o6 -27//_i,o,-i^2 - 312//_i,o,-i,o - 19//_i,o,-i,o,o 

- 1 12/9//_i,o,o + 21/2//_i,o.o6 - 54//_i,o,o,-i,o + 160//_i,o,o,o + 97/2//_i,o,o.o.o + 32//_i,o,o,i,o + 96//_i,o.i.o 
+ 8//_i,o,i,o,o - 24//_i,i^2 - 144//_i,i,_i,o - 48//_i,i,o + 96//_i,i,o,o + 48//_i,i,i,o - 17896/81//o + 493/6//o6 

- 161/18//o^2 - 97/40//o6^ - 103/2//o,-i6 - 216//o,-i^2 + 42//o,-i61n2 + 39//o,-i,-i6 - 3 12//o,- 1,-1,0 
+ 33//o,-i,-i,o,o - 1 12/9//o,-i,o - 44//o,-i,o^2 + 79//o,-i, 0,-1,0 + 102//o,-i,o,o - 99//o,-i,o,o,o - 56//o,-i,o,i,o 
+ 96//o,-i,i,o - 14//o,-i,i,o,o - 3172/27//o,o - 91/2//o,o^3 + 75//o,o^2 + 23/2//o,o,-i6 + 194//o,o,-i,-i,o 

+ 134//o,o,-i,o -251/2//o,o,-i,o,o -96//o,o,-i,i,o + 337/18//o.o,o -233/4//o,o,o6 -403/2//o,o.o.-i,o 

+ 107/2//o,o,o,o + 203/4//o,o,o,o,o + 96//o,o,o,i,o - 16//o,o,i^2 - 96//o,o, 1,-1,0 - 30//o,o,i,o + 36//o,o,i,o,o 

+ 32//o,o,i,i,o-9/2//o,i^3 + 16//o,i6 +42//o,i^2ln2-42//o,i,-i^2 + 96//o,i,-i,o- 14//o,i,-i,o,o + 32//o,i,o 

+ 38//o,],o^2 -46//o,i, 0,-1,0 + 47 //o, 1,0,0 + 63//o, 1,0,0,0 + 4O//0, 1,0, 1,0 -32//o, 1,1,0 + 26//0, 1,1, 0,0 + 44//1 ^3 - 8//i^2 

-927/20//i^2^ -24//i,_i^2 - 144//i,_i,_i,o -48//i,_i,o + 96//i,_i,o.o +48//i,_i,i,o - 16//i,o -44//i,o6 
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+ 2//i,o^2 + 75//i,o,-i^2 + 96//i,o -1,0 + 67//i,o,-i,o,o + 32//i,o,o - 1 17/2//i,o,o^2 - 42//i^o,o -1,0 - 50//i,o,o,o 
- 1 13/2i7i,o,o,o,o - 32i7i,o,i,o - 88i7i,o,i,o,o + 8i7i,i^2 +4877i,i,_i,o + 1677i,i,o - 32i7i,i,o,o - l6Hi,ixo)p(x)] 

+ eC^^|(l - 1488x5 + 3720x4 -3132x3 + 978x2 -80xi)(-2//o,-i,-iAo) + (1 - 120x5+924x4 - 1496x3 
+ 884x2 - 134xi)(47^_i,o^2) + (1 -233/2x5 +364x4 -781/2x3 + 1285/8x2 - 165/8xi)(-32i7o,o^2) 
+ (1-3312/35x5 + 1656/7x4-6708/35x3 + 1782/35x2 - 16/5xi)(-3577i,o,o^2) + (1 - 1008/1 1x5 
+ 2520/11x4 - 188x3 +582/11x2 -48/1 lxi)(22//_i,o,-i,o,o) + (1 -30096/353x5 +75240/353x4 
-60924/353x3 + 16146/353x2 - 1072/353xi)(-353/10i7i^2^) + (1 - 15792/197x5+39480/197x4 
-32076/197x3 + 8634/197x2 -640/197xi)(591/4^5) + (l -80x5 + 1352/7x4-1100/7x3 + 1018/21x2 
-8/7xi)(42i7o,i,o,o)-(l -51504/751x5 + 128760/751x4 - 105348/751x3 + 29262/751x2-2672/751x0 
(751/20//o6^) + (1 - 13680/229x5 + 34200/229x4 - 27468/229x3 + 7002/229x2 - 512/229xi) 
(229/10//-i^2^) + (l -56x5 + 140x4- 116x3 + 34x2- 4xi)(-24//o,i,o,i,o) + (l -8328/161x5 +20176/161x4 

- 16220/161x3 + 16186/483x2 - 3436/483xi )(-483/5^2^) - (1 - 2784/61x5 + 6960/61x4 - 5564/61x3 

+ 1386/61x2 - 120/61xi)(122//o,i,o,o,o) + (l - 1248/29x5 +3120/29x4 -2492/29x3 +618/29x2 -56/29xi) 

(116//o,i,o,-i,o) + (l -2384/63x5 +5960/63x4 - 14380/189x3 +410/21x2 -416/189xi)(-189/2^2^3) 

+ (1-3504/101x5 + 8760/101x4-6796/101x3 + 1434/101x2 -96/101xi)(101//o,o -1^2) + (1 -240/7x5 

+ 600/7x4-484/7x3 + 18x2 - 16/7xi)(-42//_i,o -1^2) + (1 -768/23x5 + 1920/23x4 - 1492/23x3 + 318/23x2 

-24/23xi)(-92//_i,oA-i,o) + (1-360/1 1x5 + 900/1 1x4 -700/1 1x3 + 150/1 1x2- 12/1 lxi)(-88//_i,o,i,o,o) 

+ (1-2688/85x5 + 1344/17x4- 1036/17x3 +210/17x2 -72/85xi)(34077o,o,-i,-i,o) + (l -2928/97x5 

+ 7320/97x4 - 5404/97x3 + 786/97x2 + 32/97xi)(-9777o,o,o,-i,o) + (1 - 2096/71x5 + 5240/71x4 

- 12076/213x3+798/71x2 - 176/213xi)(21377o,i^3) + (1 - 144/5x5+72x4 - 1972/35x3 +438/35x2 
-48/35xi)(210i7o,-i,-i6) + (1 -624/23x5 + 1560/23x4 - 1188/23x3+222/23x2 - 16/23xi)(138i7i,o,-i^2) 
+ (1 - 1808/67x5 +4520/67x4-10300/201x3 +630/67x2- 128/201xi)(-201//o,a-i,o,o) + (l -336/13x5 

+ 840/13x4 - 588/13x3 + 42/13x2 + 16/13xi)(26//i,o -1,0,0) + (1 - 24x5 + 6OX4 - 228/5x3 + 42/5x2 -4/5xi) 

(40//o,o,i,o,o) + (1 - 24x5 + 60x4 - 44x3 + 6x2)(32//_ 1,0,0, 1,0 - 96//o,o,- 1,1,0 - 16//o,o,i6 - 96//o,o, 1,-1,0 

+ 32//o,o,i,i,o) + (l -240/11x5 + 600/11x4 -532/11x3 + 18x2 -48/1 lxi)(-ll/2//o,o,o,o,o) + (1 - 144/7x5 

+ 360/7x4- 260/7x3 + 30/7x2)(-84//o,-i61n2 + 28//o,_i,i,o,o-84//o,i^2ln2 + 84//o,i,-i^2 + 28//o,i,-i,o,o) 

+ (1 - 20x5 + 50x4 - 1 12/3x3 + 6x2 - 2/3x1 )(-96//o,-i,o^2) + (1 - 160/9x5 + 400/9x4 - 868/27x3 + 34/9x2 

-8/27xi)(-108//o,i,i,o,o) + (l -816/49x5 +2040/49x4 - 1500/49x3 +30/7x2 -32/49xi)(-49//o,-i^3) 

+ (1 - 40/3x5 - 196/3x4 + 1528/9x3 - 84x2 - 2/3xi)(36//-i,o,o,o) + (1 - 72/7x5 + 180/7x4 - 76/7x3 - 66/7x2 

+ 20/7xi)(28i7o,o,o,i,o) + (1-164/17x5 + 372/17x4 -294/17x3 + 397/51x2 -235/51xi)(408i7o,_i^2) 

+ (1 - 28/3x5 + 2674/69x4 - 3746/69x3 + 695/23x2 - 443/69xi)(-13877o^2) + (1 - 856/121x5 + 2896/121x4 

-3372/121x3 + 1648/121x2 -42/1 lxi)(48477o,o,-i,o)-(l -80/13x5 +200/13x4-76/13x3 -86/13x2 

+ 16/1 3x1 )(39i/i,o,o,o,o) - (1 -640/173x5 + 1624/173x4 - 1444/173x3 + 826/173x2 -744/173xi)(173//o,o,o,o) 

+ (1 - 8/9x5 + 32/9x4 - 56/9x3 +59/9x2 -47/9xi)(-72i7o.o,i,o) + (1 + 644/51x5 - 318/17x4 - 154/17x3 

+ 1403/51x2 -521/51xi)(102//o,o,o) + (1 + 192/13x5 -480/13x4+420/13x3 - 150/13x2 - 8/13xi)(13//o,o6) 

+ (1 + 1416/91x5 - 336/13x4 + 488/91x3 + 684/91x2 - 282/91xi)(364//o,_i,o,o) + (1 + 792/37x5 - 1428/37x4 

+ 488/37x3 + 177/37x2 -39/37x1 )(296//i,o,o,o) + (l +720/17x5 - 1800/17x4 + 1572/17x3-558/17x2 

+ 32/17xi)(687^i,o,o,-i,o) + (1 +48x5 - 120x4 + 68x3 + 18x2 - 16xi)(-7^_i,o,o,o,o) + (1 + 1008/17x5 
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-2520/17x4 + 132x3 - 846/17x2 + 80/17xi)(34//o -ia-i,o) + (1 + 66x5 - 168x4 + 416/3x3 -467/12x2 
+ 17/12xi)(96//i,o^2) + (1 +72x5 - I8OX4 + 156x3 -54x2 + 4xi)(-16//_i, 1,0^2 - 16//^ 1,1, 0,0,0 - 16//i,_i,o^2 

- 16//i,_i,o,o,o- 16//i,o,o,i,o + 16//i,i,o^2 + 16//i,i,o,o,o) + (l +816/11x5 -2040/11x4 + 1996/11x3 
-954/11x2 + 160/llxi)(-ll/2//oA0^2) + (l + IO8X5 -270x4 + 232x3 -78x2 + 6xi)(8//o,i,o6) + (1 + 168x5 
-420x4 + 364x3 - 126x2 + 12xi)(87/_i,o^3) - (1 + 2064/5x5 - 1032x4 + 4404/5x3 - 1446/5x2 + 128/5xi) 
(107/0-1,0,0,0) + (1 + 528x5 - 1320x4 + 1 124x3 - 366x2 + 32xi)(-//_i,o,o^2) + (1 + 5392x5 - 13960x4 + 1 1796x3 
- 3162x2 + 68xi)(-//o^3) + (1 -48x4+4320/37x3 - 3120/37x2 +466/37xi)(-74i/_i^2) + (1 -4184/235x4 

+ 9152/235x3 - 1078/47x2 +298/235xi)(235^2) - (1 - 186/13x4 + 372/13x3 - 177/13x2 -9/13xi) 
(208//_i _i,o,o) - (1 - 63/5x4 + 126/5x3 - 21/2x2 - 21/10xi)(1607/_i,o,_i,o) + (1 + 24/13x4 - 240/13x3 
+ 308/13x2 -92/13xi)(52i/i,o,o) + (1 + 336/137x4 + 960/137x3 -2080/137x2 + 674/137xi)(2747/_i,o,o) 
+ (1 + 208/37x4 - 608/37x3 + 604/37x2 - 204/37xi)(148//o,i,o) + (1 + 8x4 - 16x3 + 143/15x2 - 23/15xi) 
(240//i,i,o,o) + (l +204/25x4 -408/25x3 + 372/25x2- 168/25xi)(100//-i ^3) + (1 + 1184/139x4 
-4000/139x3+4280/139x2 - 1354/139xi)(-278//o,_i,o) + (l +2152/213x4 - 1872/71x3+5170/213x2 

- 1714/213xi)(213//o,o) + (l + 12x4 -24x3 + 32/3x2 + 4/3xi)(144//_i^2ln2 -48//-i, 1,0,0 + 144//i^2ln2 

- 144//i,_i^2-48//i,_i,o,o) + (l + 14x4-28x3 + 167/9x2 -41/9xi)(-432//_i,_i^2) - (1 +7296/395x4 

- 17616/395x3 +2592/79x2 -462/79xi)(790/3^3) -(1 + 1572/79x4-3144/79x3 + 1746/79x2 - 174/79xi) 
(316//i^3) + (1 +21x4-42x3 + 85/4x2 - l/4xi)(64//i,o,i,o) - (1 + 156/5x4 - 312/5x3 + 65/2x2 - 13/lOxi) 
(160//i,o,-i,o) + (1 + 32x4 - 64x3 + 767/21x2 - 95/21xi)(-336//o,-i,-i,o) + (1 + 48x4 - 96x3 

+ 54x2 - 6xi)(487/o -1,1,0 + SHo,i^2 + 4SHox-i,o - 167/o,i,i,o) - (1 + 72x4 - 144x3 + 132x2 - 60xi) 
(24^2ln2) + (1 + 96x4 - 192x3 + 100x2 - 4xi)(-167/_i,o,i,o) + (1 - 8x3 + 12x2 - 6x1) 
(44//i,o) - (1 +240/187x3 -360/187x2 -254/187xi)(187/4i/o) + (1 + 1344/55x3-2016/55x2 
+ 562/55x1 )(-2207/_i,_i,o) + (1 + 568/9x3 - 284/3x2 + 266/9xi)(367/_i,o) + (1 + 6x2 - 6xi)(81n'*2 
+ 192Li4(l/2) + 96^2ln^2) + (4 - 20x2 + 20xi) + (x5 - 5/2x4 + 2x3 - l/2x2)(768//i,o^3 + 1536//i,o,i,o,o) 
+ (x5 - 5/2x4 + 13/6x3 - 3/4x2 + l/12xi)(-1536//o,-i,o,i,o) + {x3 - 3/2x2 + 1/2x0(1 152//_i,i,o 
+ 192//i^2 + 1152//i,_i,o -384//i,i,o) + (-369/4^5 +400^3 + 52^2-29/2^2^3 - 175^2^ + I6//-16 
-240//_i^2 -731/10//_i^2^ - 16//-i,-i^2 - 32//_i,_i,_i,o - 1248//_i,_i,o + 16//_i._i.o,o + 360//_i,o 
-56//^i,o^3 - 44//_i,o^2 + 54//_i,o,-i6 + 16//-i,o,-i,o + 54//_ 1,0,- 1,0,0 + 816//_i,o,o - 65//_i,o,o6 
+ 4//-i,o,o,-i,o -52//_i,o,o,o - 81//- 1,0,0,0,0 - 32//_ 1,0,0, 1,0 -5 6//_ 1,0, 1,0,0 + 384//_i,i,o - 16//_i, 1,0^2 

- 16//_i,i,o,o,o + 73//0 + 169/3//o^3 + 103//o^2 + 1633/20//o^2^ + 1 1 l//o,-i^3 + 320//o,-i^2 - 84//o,_i^2ln2 

- 78//o,-i,-i6 + I6//0.-1.-1.0 - 98//o,-i,-i,o,o + 816//o,-i,o + 64//o._i.o^2 - 190//o.-i.o.-i.o + 6O//0.-1.0.0 

+ 214//o,_i,o,o,o + 128^0,-1,0,1,0 + 28^o,_i, 1,0,0 - 244//o,o + 37^o,o^3 - 192//o,o^2 + 37//o,o,-i^2 - 76^o,o,-i,-i,o 
+ 20//o,o,-i,o + 55//o,o,-i,o,o + 96//o,o,-i,i,o - 313//o,o,o + 309/2//o,o,o^2 + 351//o,o,o,-i,o - 197//o,o,o,o 

- 123/2//o,o,o,o,o - 180//o,o,o,i,o + 16//o,o,i^2 + 96//o,o, 1,-1,0 - 8//o,o,i,o + 40//o,o,i,o,o - 32//o,o,i,i,o + 69//o,i^3 

- 84//o,i^2ln2 + 84//o,i,_i^2 + 28/?o,i,-i,o,o - 256/?o,i,o - 40/?o,i,o6 + 2 12/?o,i, 0,-1,0 - 214/?o,i,o,o 

- 170//0.1.0.0.0 - I2O//0.1.0.1.0 - 12//o,i,i,o,o + 64//i^2 + 227/2//i^2^ + 384//i,_i,o - 16//i,_i,o6 - 16//1,- 1,0,0,0 

- 128//i,o + 96//i,o^3 - 128/?i,o6- 150/^1,0,-1^2 - 1 34//i, 0,-1, 0,0 -256//i, 0,0 + 141//i,o,o6 + 132//i,o.o,-i,o 

- 128//i,o,o,o + 105//i,o,o,o,o - 16//i,o,o,i,o + 192//i,o,i,o,o - 128//i,i,o + 16//i,i,o^2 + 16//i,i,o,o,o)K^) 
+ (222^5-32^3-92^2^3 + 192/5^^ + ^4//^^^^ -488/5//-i62 _256//_j _j _ 64//_j ^,^2 

- 128//_i,_i,o,-i,o + 256//_i,_i,o,o + 192//_i,_i,o,o,o + 128//_i,o^2 - 32//_i,o,-i^2 - 64//_i,o,-i,-i,o 
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+ 160//- 1,0 -1,0,0 - 256//_i,o,o - 96//-i,o,o^2 + 160//- 1,0,0,-1,0 - 64//_i,o,o,o - 272//_ 1,0,0,0,0 - 64//_i,o,o,i,o 
+ 128//_i,o,i,o-64//_i,i^3 -64//_i,i,o^2 + 128//_i,i,o,-i,o - 128//_ 1,1, 0,0,0 

- 128//-i,i,o,i,o - 16//o^3 -96//o6 - V02/5//o^2' - 16//o,-i^3 + 64//o,-i^2 - 16//o,-i,-i^2 

- 32//o,_i,_i,_i,o - 256//o,-i,-i,o + 144//o,-i,-i,o,o - 64//o,-i,o + 48//o,-i,o^2 + 144//o,-i,o,-i,o 

+ 96//o,-i,o,o - 120//o,-i,o,o,o - 64//o,-i,o,i,o + 192//o,-i,i,o + 192//o,o - 296/3//o,o6 + 128//o,o^2 + 152//o,o,-i^2 
+ 208//o,o,-i,-i,o + 224//o,o,_i,o - 136//o,o,-i,o,o-96//o,o,-i,i,o + 128//o,o,o - 176//o,o,o^2 -408//o,o,o,-i,o 
+ 80//o,o,o,o + 1 16//o,o,o,o,o + 160//o,o,o,i,o - 16//o,o,i^2 - 96//o,o,i,_i,o - 192//o,o,i,o - 80//o,o,i,o,o + 32//o,o,i,i,o 
+ 64//o,i^3 + 32//o,i^2 + 192//o,i,-i,o - 64//o,i,o - 32//o,i,o^2 - 128//o,i, 0,-1,0 ~ 128//o,i,o,o + 32//o,i,o,o,o 
+ 128//o,i,o,i,o - 64//o,i,i,o + 296/5//i^2^ - 64//i,_i^3 - 64//i,_i,o^2 + 128//i,_i, 0,-1,0 - 128//i,_i,o,o,o 

- 128//i,_i,o,i,o + 16//i,o^3 - 32//i,o,-i^2 + 320//i, 0,-1, -1,0 - 256//i,o,-i,o,o ~ 192//i,o,-i,i,o - 32//i,o,o^2 

- 384//i,o,o,-i,o + 192//i,o,o,o,o + 256//i,o,o,i,o - 32//i,o,i^2 - 192//i,o,i,-i,o + 128//i,o,i,o,o 

+ 64//i,o,i,i,o + 64//i,i^3 + 64//i,i,o^2 - 128//i,i, 0,-1,0 + 128//i,i,o,o,o + nSHixoxo)p(xf) 

CT;Cpr^|5204/8 1 + (1 - 2x1 )(8//-i6 + 16//-i,-i,o - 8/9//_i,o - 8//-i,o,o - 203/27//o + 2//o^2 - 8//o,-i,o 

+ 4/9//o,o + 4//„,„,o) + (l + l/39xi)(104/3^2) + (l +6x1X8/3^3) + (-448/9^3 - 1136/27^2 -96/5^^ 
+ 32//_i^3 + 448/9//_i^2 - 32//_i,_i^2 - 64//_i,_i,_i,o + 896/9//_i,_i,o + 32//_i,_i,„,o - 2272/27//^i,o 

- 8//_i,o^2 + 32//_i,o,-i,o - 448/9//_i,o,o - 16//_i,o,o,o + 5204/81//o - 16/3//o^3 + 92/9//o^2 + 16//o,-i^2 
+ 32//o,-i,-i,o - 448/9//o,-i.o - I6//0.-1.0.0 + 1 136/27//o,o + 4//o,o^2 - 16//o,o,-i,o + 224/9//o.o.o 

+ 8//o,OAo)pW}+^5(e^), 

^2^^' = (22) 
C^{-2^(x)//o} 

+ 6Cp|^(x)(2^2 + 4//-i,o - 8//0 - 2//o,o)} 

+ e2c^{^(x)(4^3 + 8^2 - 4//_i6 - 8//_i,_i,o + 16//_i,o + 4//-i,o,o - 16//0 + //q^ + 4//o,-i,o 
-8//o,o-2//o.o.o)} + ae^), 



_ ^2)//o + 4^(x)//o + Sq(xfHo,o] (23) 

C^r^ |96(X5 - 5/2x4 + 7/4x3 - 1/8x2 - l/8xi)(//o^2 + //o,o,o) - 544^2(.^4 - 2x3 + 33/34x2 + l/34xi) 

- 352/3//o,o(x4 - 2x3 + 23/22x2 - l/22xi) - 496/3//o(x3 - 3/2x2 + l/2xi) + 272/3(xi - X2) 

- (272/9//0 + 6//o^2 + 6//o,o.o)^w} 

C^|l032//oAoU5 - 136/43x4 + 603/172x3 - 1555/1032x2 + 169/1032x0 + 1032//o6fe - 124/43x4 
+ 507/172x3 - 1267/1032x2 + 169/1032x0 - 12//o,o.o,ofe - 5/2x4 + 13/12x3 + 7/8x2 - 1 1/24x0 

- 192//o,o,i,oU5 - 5/2x4 + 7/4x3 - 1/8x2 - 1/8x0 + 360//o,_i^2(x5 - 5/2x4 + 1 13/60x3 - 13/40x2 - 7/120x0 
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- 1296//o^3fe - 5/2x4 + 635/324x3 - 95/216x2 - 13/648xi ) + (x5 - 5/2x4 + 55/28x3 - 25/56x2 

- 1/56x1 )1344(//i,o6 + //i ,0,0,0) + (X5 - 5/2x4 + 547/276x3 - 87/184x2 - 5/552xi )(276//o,o^2) 

+ (X5 - 5/2x4 + 1 123/564x3 - 183/376x2 - 5/1 128xi )(1 128//o,_i,o,o) + (xj - 5/2x4 + 2419/1 172x3 

- 1397/2344x2 +75/2344xi)(1758/5^2^)-96//o,i,o,ofe -5/2x4 + 35/12x3 - 15/8x2 + ll/24xi) 
-48//o,o-i,oU5 -5/2x4 + 49/12x3 -29/8x2 + 25/24xi) + (x4- 121/56x3 +485/336x2 -95/336xi)(672^2) 

- 192i7o,i,o(^4 - 2x3 + X2) - 1296^3(x4 - 2x3 + 86/81x2 - 5/81xi) + (x4 - 2x3 + 52/47x2 - 5/47xi) 

(1 128//_i,o,o) + (X4 - 2x3 + 4/3x2 - l/3xi)(360i7_i^2) -48(77o.-i.o + 2i/i,o.o)(x4 - 2x3 + 5/3x2 - 2/3xi) 
-96Ho,o(xA - 13/8x3 + 13/48x2 + 17/48xi) - 12(277_i,o + 3Ho + 8/?i,o)(x3 - 3/2x2 + l/2xi) 
- 12(1 + 8^2ln2)(x2 -xi) + (1/3^2+45/8^2^ + 2/3i?_i,o -254/9//o + 13i/o^3 + 17/2//o6 -21/277o,-i^2 
-5/2^0,-1,0,0 - 1/3^0,0 - 5/4^0,0^2 - 25i?o,o,-i,o + 17/277o,o,o + 1 l/4^o,o,o,o + 12i7o,o,i,o - 22^o,i,o,o 

- 12i7i,o^2 - 12//i,o,o,o)?W + (12^2 - 18i/-i6 - 2i7_i,o,o + IHo^i - 4i/o,-i,o + 12i/o,o + 3Ho,o,o)q(xf 

+ (69/10^2^ + 12//06 - 18i/o,-i^2 - 277o,-i.o.o + 1Hofi^2 - 4i/o,o,-i,o + 12i/o,o,o + 3i/o,o,o,o)?W^) 



C/\^(X5 - 18/5x4 + 251/60x3 - 169/120x2 -7/40xi)(240//o,o,o) + (xj - 13/5x4 + 131/60x3 

- 1/120x2 -23/40xi)(240//o6) + fe - 5/2x4 + 13/12x3 +7/8x2- ll/24xi)(24//o,(),(),()) + (x5 -5/2x4 
+ 7/4x3 - 1/8x2 - 1/8x1 )(384//o,o,i,o) - 1536//o,o,-i,o(x5 -5/2x4 + 89/48x3 -9/32x2 -7/96x0 

- 14407/o,i,o,o(.X5 -5/2x4 + 67/36x3 -7/24x2 -5/72xi)- 7207^0,-1^2(^5 -5/2x4 + 113/60x3 - 13/40x2 

- 7/120x0 + (X5 - 5/2x4 + 2171/1 140x3 - 271/760x2 - 109/2280x0(684^2^) + (^5 - 5/2x4 + 71/36x3 

- 11/24x2 - l/72xi)(1008//o,_i,o,o) + (X5 -5/2x4 + 263/132x3 -43/88x2 - 1 /264xi)(26477o,o^2) 

+ (X5 - 5/2x4 + 25/12x3 - 5/8x2 + l/24xi)(-67277o6 + 57677i,o6 + 57677i,o,o,o) - 76877o,oU4 - 21/8x3 
+ 37/16x2 - 11/16x0 + (^4 - 11/5x3 + 19/15x2 - 1/15x0(960^2) - 1440//i,o,oU4 -2x3 +44/45x2 
+ 1 /45xi) + (X4 - 2x3 + 64/63x2 - l/63xi)(1008i7_i,o,o) - 15367^0 _i,o(x4 - 2x3 + 49/48x2 - 1 /48xi) 
-672^3(X4 - 2x3 +43/42x2 - 1/42x0 + (x4 - 2x3 + 13/12x2 - l/12xi)(384//o,i,o) 

-720//_i^2(^4 - 2x3+4/3x2 - 1/3x0 + (X3 - 3/2x2 + l/2xi)(-768//_i,o + 192//i,o) + (x2 - xi)(192^2ln2) 
+ (-406 + 36//-16 + 4//_i,o,o - 38//o6 - 8//o,-i,o + 8//o,o - 14//o,o,o - 16//o,i,o - l6Hisio)qixf 
+ (-34^ - 327/20^2 _ ^ + 31^^, _ ^4Ho^3 + 43//o6 + 21//o,-i6 - lHo,-i,o,o + 10//o,o - I/2//0.06 
+ 56//o,o,-i,o + 43//o,o,o - 1 1 /2//o,o,o,o - 24//o,o,i,o + 50//o,i,o,o - 16//i,o + 12//i,o6 + 12//i,o,o,o)^(x) 
+ (-69/5^2 + 36//o,-i6 + 4//o,_i,o,o - 14//o,o6 + 8//o,o,-i,o - 6Ho,o,o,o)q(x?] 



eCpTj^l^ixs - 31/6x4 + 1031/132x3 - 1057/264x2 + 95/264xi)(1767/o,o,o) + (X5 - 23/6x4 + 679/132x3 

- 689/264x2 + 79/264xi)(176//o6) + fe - 5/2x4 + 7/4x3 - 1/8x2 - l/8xi)(-48/562 + 192//_i,o6 

+ 192//_i,o,o,o - 192//06 - 576//o,-i6 - 192//o,-i,o,o + 96//o,o6 - 192//o,o,_i,o + 288//o,o,o,o + 384//o,i,o,o) 

+ (X4 - 605/287x3 + 49/41x2 - 25/287xi)(-4592/36) + (X4 - 2x3 + 55/59x2 + 4/59xi)(-3776/36) 

+ (X4 - 2x3 + 19/20x2 + l/20xi)(256061n2) + (X4 - 2x3 + 23/22x2 - l/22xi)(704//_i6 + 704/37/_i,o,o 

+ 704/3//o,-i,o - 1408/3//i,o,o) + (X4 - 205/149x3 + 34/ 149x2 + 22/ 149xi )(-2384/9//o,o) 

+ (X3 - 3/2x2 + l/2xi)(992/3//_i,o - 4040/9//o) + (-1528/9x2 + 1528/9x0 + (272/9^ + 3/5^^ 

+ 544/9//_i,o - 12//_i,o6 - I2//-1.0.0.0 - 1 m/llHo + UHo^s + 59/3//o6 + 36//o,-i^2 



^iCpTj (-8/3^ - 16/3//_i,o + IOO/9//0 + 8/3//o,o)^(x) 
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+ 12//o,-i,o,o - 272/9i/o,o - 6//o,o^2 + 12i/o,o,-i,o + 55/3//o,o,o - 18i/o.o.o.o - 24//o,i,o,o)^W| 

+ eCpC^|(;c5 -916/247x4 +4663/988x3 -4469/1976x2 +495/1976xi)(-3952//o,o,-i,o) 
+ (X5 - 521/146x4 + 1263/292x3 - 1 125/584x2 + 99/584xi)(2336//o,o,i,o) + - 144/41x4 
+ 701/164x3 -695/328x2 + 117/328xi)(-1968//o,-i6) + (^5 - 1592/491x4 + 7107/1964x3 
-5857/3928x2 +45 l/3928xi)(-392877o,o6) + Us -2329/730x4 +4911/1460x3 -3393/2920x2 
-33/2920xi)(1460i7o,o,o) + Us -49/16x4 + 1209/368x3 - 1021/736x2 + 121/736xi)(147277o,o,o,o) 
+ (X5 - 1217/407x4+4959/1628x3-3637/3256x2 + 199/3256xi)(-651277o,-iAo) + Us -2173/730x4 
+ 4647/1460x3-4137/2920x2 + 123 /584xi)(146077o6) + Us- 1295/458x4 + 7429/2748x3 

- 1705/1832x2 + 301/5496xi)(5496//o^3)- Us - 1367/530x4+479/212x3-309/424x2 + 103/2120xi) 
(8480//i,o^2) - Us - 285 1 / 11 17x4 + 49409/22340x3 - 3093 1 /44680X2 + 1473/44680x1 )(4468^2^) 

+ (X5 - 5/2x4 + 7/20x3 + 79/40x2 - 33/40xi)(-480//(),(),-i,-i,()) + (xj - 5/2x4 + 161/156x3 + 99/104x2 

- 151/312xi)(156^2^3) + Us -5/2x4 + 13/12x3 +7/8x2- ll/24xi)(-24//_i,o,o,o,o) + Us -5/2x4 

+ 457/348x3 + 123/232x2 -239/696xi)(-696//o,i^3) + Us -5/2x4 + 157/108x3 +23/72x2 -59/216xi) 
(648//o,o,o -i,o) + Us - 5/2x4 + 139/84x3 + 1/56x2 - 29/168xi)(-336//i,o -1,0,0) + Us - 5/2x4 

+ 1643/948x3 -63/632x2 -253/ 1896x1 )(1896//o,o,-i,o,o) + Us -5/2x4 + 7/4x3 - 1/8x2-1/8x1) 
(-384//_i,o,o,i,o + 864//o,-i61n2 - 288//o,-i,i,o,o + 1 152//o,o,-i,i,o + 192i/o,o,i6 + 1 152//o,o, 1,-1,0 - 384i7o,o,i,i,o 

+ 864i7o,i^2ln2 - 864i7o,i -1^2 -28877o,i -1,0,0) - Us -5/2x4 + 281/156x3 - 21/104x2 - 31/312xi) 
(1872//i,o,-i^2) - 4807?o,o,i,o,oUs - 5/2x4 + 109/60x3 - 9/40x2 - 1 1/120x0 - 960//o,o,o,i,oUs - 5/2x4 

+ 11/6x3 - 1/4x2 - l/12xi) + (X5 -5/2x4 + 155/84x3 - 15/56x2 - 13/168xi)(-3024i7o,-i,-i6) 

+ (X5 - 5/2x4 + 1 13/60x3 - 13/40x2 - 7/120xi)(720i7-i,o,-i^2) + Us - 5/2x4 + 61 1/324x3 - 71/216x2 

- 37/648x1 )(2592//o,i,i,o,o) + Us - 5/2x4 + 205/108x3 - 25/72x2 - 1 l/216xi)(2592//o,-i,o^2) 

+ (X5 - 5/2x4 + 1669/876x3 - 209/584x2 - 83/1752xi)(1752//i,o,o,o,o) + Us - 5/2x4 + 389/204x3 - 49/136x2 

- 19/408xi)(1020//o,o,o^2) + Us -5/2x4 + 1961/1020x3 -261/680x2 -79/2040xi)(2040//o,-i^3) 
+ (X5 - 5/2x4 + 539/276x3 - 79/184x2 - 13/552xi )(-5520//o,-i,o,-i,o) + Us - 5/2x4 + 305/156x3 

- 45/104x2 - 7/312xi)(2496//o,-i,o,i,o) - Us - 5/2x4 + 635/324x3 - 95/216x2 - 13/648xi )(2592//_i,o^3) 

+ (X5 - 5/2x4 + 55/28x3 - 25/56x2 - l/56xi)(2688//-i,i,o^2 + 2688//_i, 1,0,0,0 - 2016//o,i,o6 + 2688//i,-i,o^2 
+ 26887/1,-1,0,0,0 + 26887^1,0,0,1,0 - 2688//1, 1,0^2 - 26887^ij,o,o,o) + Us - 5/2x4 + 2051/1044x3 

- 3 1 1 /696x2 - 37/2088xi)(-8352//i,o,o,-i,o) + Us - 5/2x4 + 3797/1932x3 - 577/1288x2 - 67/3864x0 
(77287^0,-1,0,0,0) - (X5 - 5/2x4 + 1 1 1 1 /564x3 - 171/376x2 - 17/1 128x0(1 12877o,o6) - Us - 5/2x4 

+ 71/36x3 - 11/24x2 - l/72xi)(50407?o,-i,-i,o,o) + (xs -5/2x4 + 547/276x3 -87/184x2 -5/552x0 
(55277_i,o,o6) + Us - 5/2x4 + 1 19/60x3 - 19/40x2 - l/120xi)(-96077o,i,o,i,o) + (^s - 5/2x4 + 1 123/564x3 

- 183/376x2 -5/1128xi)(225677_i,o,-i,o,o) + (^s -5/2x4 + 313/156x3 -53/104x2 + 1/312x0 
(-12487/1,0^3 - 24967/1,0,1, 0,0) + (^s - 5/2x4 + 5001/2492x3 - 2543/4984x2 + 17/4984xi )(1 1214^5) 

- 32407/i,o,o^2(xs - 5/2x4 + 3259/1620x3 - 559/1080x2 + 19/3240x0 + (xj - 5/2x4 + 9195/4564x3 
-4765/9128x2 + 67/9128xi)(-13692/57/i^2^) + (x5 -5/2x4 + 3347/1652x3 -1781/3304x2 

+ 43/3304xi)(-7434/57/o^2^) - (^s - 5/2x4 + 1501/732x3 - 281/488x2 + 37/1464xi)(29287/o,i,o,o,o) 
+ (X5 - 5/2x4 + 617/300x3 - 1 17/200x2 + 17/600xi )(24007/o,i,o,_i,o) + (xj - 5/2x4 + 2419/1 172x3 

- 1397/2344x2 + 75/2344xi)(3516/57/_i^2^) + (^s -5/2x4 +43/20x3 -29/40x2 + 3/40x0 
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(60//o,o,o,o,o) + - 5/2x4 + 267/ 1 16x3 - 221/232x2 + 35/232xi)(6967/o,o -1^2) 
+ (X5 -5/2x4 + 35/12x3 - 15/8x2 + ll/24xi)(-192//_i,(),i,o,o) 

-96//-i,o,o,-i,ofe -5/2x4 + 49/12x3 -29/8x2 +25/24x0- 8480//i,o,o,ofe - 131/53x4 + 2167/1060x3 

- 1261/2120x2 +47/2120x0 + (X5 - 163/86x4 + 167/172x3-67/1032x2 - Il/1032x0(2064//_i,o6) 
+ (X5 - 131/86x4 - 83/172x3 +453/344x2 - 107/344xO(688//o,iAo) + + 51/86x4 - 689/172x3 

+ 911/344x2 -81/344xO(20647?_i,o,o.o) + U4 -428/121x3 +999/242x2 -35/22xO(4847/o,o) 
+ (X4 - 122/49x3 + 97/49x2 - 24/49xO(784//i,o,o) + (^4 - 9/4x3 + 391/216x2 - 121/216x0(864^3) 
+ (X4 - 574/257x3 + 885/514x2 - 251/514x0(1028^2) + (xa - 2x3 + 3/5x2 + 2/5xO(-48077o,-i,-i,o) 
+ (X4 - 2x3 +21/29x2 + 8/29xO(-696i7i^3) + (x4 - 2x3 + 7/9x2 + 2/9xO(86477_i^2ln2 

- 28877_i,i,o,o + 8647^1 ^2ln2 - 864//i,_i^2 - 288i7i,_i,o,o) + (x4 - 2x3 + X2)(l 152i/o,-i,i,o 
+ 192//o,i6 + 1 152//o,i,-i,o - 384//o,i,i,o) + (X4 - 2x3 + 85/81x2 - 4/81xO(2592/?i4,o,o) 

+ (X4 - 2x3 + 85/78x2 - 7/78xO(2496//_i,o,i,o) - (^4 - 2x3 + 346/315x2 - 31/315xO(5040//_i,_i,o,o) 
+ (X4 - 2x3 + 76/69x2 - 7/69xO(-5520//_i,o,_i,()) + (X4 - 2x3 + 17/15x2 - 2/15xi )(-960//i, 0,1,0) 
+ (X4 - 2x3 + 99/85x2 - 14/85x1 )(2040//_i ^3) + (m - 2x3 + 88/75x2 - 13/75xi )(2400//i,o,-i,o) 
+ (X4 - 2x3 + 74/63x2 - 1 1/63x1 )(-3024//_i,-i^2) + (m - 2x3 + 5/3x2 - 2/3xi)(-864^2ln2) 
+ (X4- 129/74x3 + 125/148x2 - 15/148xi)(1776//-i^2) + (x4- 317/184x3 +263/368x2 
+ 3/368xi)(1472//o -1,0) + (•^4 - 145/98x3 + 33/196x2 +61/196xi)(7847^_i,o,o) 

- (X4 - 22/17x3 + 5/17xO(5447/o,i,o) + (X3 - 3/2x2 + 1 /2xO(-2407/_i _i,o + 83277_i,o + 576//_i,i,o -imHo 
+ 96771^2 + 576i7i,_i,o - 65677i,o - 19277i,i,o) + (xi - xi)i-lS + 161n^ 2 + 384Li4(l /2) + 192^2ln^ 2) 

+ (153/8^5 + 2^3 + 295/9^2 - 15 1/4^2^3 - 297/20^2^ - 277-i^2 + 45/477_i^2^ - 4H_i,-i,o + 374/9i7_i,o 
+ 26i7_i,o^3 + 17i?-i,o^2 - 21i/_i,o,-i^2 - 5i/- 1,0,- 1,0,0 + 2i/-i,o,o 

- 5/2/?_i,o,o^2 - 50/?_i,o,o,-i,o + 17i?-i,o,o,o + 11 /2i?- 1,0,0,0,0 + 24i7_i,o,o,i,o - 4477_ 1,0, 1,0,0 - 24i7_i4,o^2 
-24//_i,i,o,o,o-5446/27//o + 13/2//o6 -21/4//o^2 -387/40//o6^-79/2//o,-i^3 - 15//o,-i6 

- 54//o,-i^2ln2 + 1 17//o,-i,-i^2 + 35//o,- 1,-1, 0,0 + 2Ho,-u) - 66//o,-i,o^2 + 65//o,_i, 0,-1,0 + 25//o,-i,o,o 

- 67//o,-i,o,o,o - 28//o,-i,o,i,o + 18^o,-i,i,o,o - 133/9//o,o + 17/2//o,o^3 - 91/2//o,o^2 + 105/2//o,o,-i^2 

+ 198//o,o,-i -1,0 - 151//o,o,-i,o -253/2//o,o,-i,o,o -72//o,o,-i,i,o-9/4//o,o,o -95/4//o,o,o^2 - 177/2//o,o,o,-i,o 

- 3//o,o,o,() + 9/4//o,o,o,o,o + 40//o,o,o,i,o - 12//o,o,i^2 - 72//o,o,i,-i,o + 54//o,o,i,o + 22//o,o, 1,0,0 + 24//o,o,i,i,o 
+ 239/2//o,i^3 - 54//o,i^2ln2 + 54//o,i,-i6 + 18//o,i,-i,o,o + 18^0,1,0^2 + 34//o,i,o,-i,o - 123//o,i,o,o 

- 37//o,i,o,o,o + 4//o,i,o,i,o - 74//o,i, 1,0,0 - 201/20//i^2^ -24//i,_i,o^2 - 24//i,_i,o,o,o + 12//i,o - 2//i,o^3 
-46//i,o^2 + 93//i,o,_i^2 + 29//i.o,-i.o,o - 19/2//i,o.o6 + 74//i,o,o,-i,o 

- 46//i,o,o,o - 83/27/1,0,0,0,0 - 24//i,o,o,i,o - 4i/i,o,i,o,o + 24//i,i,o^2 + 24i7i,i,o,o,o)?(Jc) 

+ (39/2^5 - 45^26 + 93/5^2^ + 69/5i/_i^2^ + 24//_i,o^2 - 3677_i,o,-i6 - 4i?- 1,0,- 1,0,0 + 14i/-i,o,o6 

- 8//_i,o,o,-i,o + 2477_i,o,o,o + 677_i,o,o,o,o - 14i/o6 +4477o6 - 307/10//o^2^ - 34i7o,-i^3 - 13277o,-i6 
+ 24i/o,-i^2ln2 + 8477o,-i,-i^2 + 12i/o,- 1,-1, 0,0 - 32i7o,-i,o^2 + 20i7o,-i, 0,-1,0 - 12^o,-i,o,o - 2077o,- 1,0,0,0 

- 8/?o,-i,i,o,o + 6i/o,o^3 + 14i?o,o^2 - 6i/o,o,-i^2 + 247^0,0,-1,-1,0 - 10477o,o,-i,o 

- 18//o,o,- 1,0,0 + 44//o,o,o - 3//o,o,o^2 - 18i?o,o,o,-i,o + 32//o,o,o,o + ^0,0,0,0,0 + 24//o,o,i,o + 2//o,i^3 

+ 24//o,i^2ln2 - 24//o,i,-i^2 - 8//o,i,-i,o,o + 8//o,i,o^2 - 8//o,i,o,-i,o - 28//o,i,o,o + 4//o,i,o,o,o - 8//o,i, 1,0,0 

- 69/5//i^2^ - 24//i,o6 + 36//1.0.-1 ^2 + 4//i,o,-i,o,o - 14//i,o,o^2 + 8//i,o,o,-i,o - 24//i,o,o,o - 6//i,o,o,o,o)^(x)^ 
+ (-12^3 + 68^2 + 24^2ln2 - 169/10^2^ - 347^_i^3 - 1927^_i^2 + 247^_i^2ln2 + 847^_i,_i^2 
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+ 12//_i,_i,o,o - 32//_i,o6 + 20//-i,o,-i,o - 24//_i,o,o - 20//_i,o,o,o - 8//-i,i,(),o + 6//o^3 + - 42//o,-i^2 

+ 24//o,-i,-i,() - 120//()^_ij) - 22//0 _i_o_o + 32//()j) + 1 l//o,o^2 - 26//o,o,-i,o + 66//()_o_o + V//o,o,o,o + 24//o_i_o 

+ 2//i^3 + 24//i^2ln2 - 24//i,_i^2 - 8//i,_i,o,o + 8//i,o^2 - 8//i,o,-i,o - 36//i,o,o + 4//i,o,o,o - 8//i,i,o,o)^W^} 



+ €C/!^(X5 - 77/10x4 + 743/60x3 - 817/120x2 + 9/8xi)(48077_i,o^2) 

+ (X5 - 191/46x4+23/4x3 -4345/1288x2 + 999/1288xi)(-128877o6) + (^5 -4x4 + 29/4x3 -75/8x2 
+ 41/8xi)(-64//o,o,i,o) + U5 -362/107x4 + 1669/428x3 - 1511/856x2 + 21 3/856xi)(34247?oA-i,o) 
+ (X5 - 728/233x4 + 3075/932x3 - 2225/1864x2 + 35/1864xi)(-372877o,o^2) 

+ (X5 - 1745/674x4 + 2847/1348x3 - 1265/2696x2 - 145/2696xi)(5392i7o6) + - 28/1 1x4 + 89/44x3 
- 367/792x2 - l/72xi)(-6336//i,o^2) + (X5 -203/80x4 + 179/80x3 - 143/160x2 + 31/160xi)(-6407/o,o,o,o) 

- (X5 - 5/2x4 + 49/60x3 + 51/40x2 - 71/120xi)(240//i,o,o,o,o) 

+ (X5 - 5/2x4 + 35/36x3 + 25/24x2 - 37/72xi)(288//(),(),(),i,()) + (X5 - 5/2x4 + 13/12x3 + 7/8x2 

- 11/24x1 )(48//-i,o,o,o,o) + (^5 -5/2x4 + 139/84x3 + 1/56x2 -29/168xi)(6727/i,o -1,0,0) 

- (X5 - 5/2x4 + 341/204x3 - 1/136x2 - 67/408xi)(816//o -16) 

+ (X5 -5/2x4 + 103/60x3-3/40x2 - 17/120xi )(-1920//o,i, 1,0,0) + fe - 5/2x4 + 7/4x3 - 1/8x2 - l/8xi) 
(768//_i,o,o,i,o - 1728//o,-i61n2 - 1920//o,-i,o6 + 576//o,-i,i,o,o-2304//o,o,-i,i,o -2928//(),o,o,-i,o 

- 384//o,o,i^2 - 2304//o,o,i,-i,o + 768//o,o,i,i,o - 1728//o,i^2ln2 + 1728//o,i,-i^2 + 576//o,i,-i,o,o) 

+ (X5 - 5/2x4 + 28 1 / 156x3 - 2 1 / 1 04x2 - 3 1 /3 1 2xi )(3744//i,o,- 1 ^2) + - 5 /2x4 + 2453 / 1 356x3 - 1 93 /904x2 

-259/27 12x1 )(-542477o,o,-i,o,o) + (X5 -5/2x4 + 109/60x3-9/40x2 - ll/120xi)(96077o,o,i,o,o) 

+ (X5 - 5/2x4 + 2887/1572x3 - 267/1048x2 - 257/3144xi)(628877o,i^3) + - 5/2x4 + 103/56x3 

-29/112x2 -9/1 12xi)(10752i?o,o,-i,-i,o) + (X5 -5/2x4 + 155/84x3 - 15/56x2 - 13/168xi)(604877o,-i,-i6) 

+ (X5 - 5/2x4 + 89/48x3 - 9/32x2 - 7/96xi)(-307277_i,o,o,-i,o) + - 5/2x4 + 543/292x3 - 169/584x2 

-41/584xi)(3504i?o.o -1^2) + (^5 - 5/2x4 + 67/36x3 - 7/24x2 - 5/72xi)(-2880//_i,o,i,o,o) 

+ (X5 - 5/2x4 + 1 13/60x3 - 13/40x2 - 7/120xi)(-1440//_i,o,-i6) + - 5/2x4 + 3401/1788x3 - 421/1 192x2 

- 175/3576x1 )(-3576^2^3) + - 5/2x4 + 2171/1 140x3 - 271/760x2 - 109/2280xi )(1368//_i ^2^) 

+ (X5 - 5/2x4 + 23/12x3 - 3/8x2 - l/24xi)(4992//o,i,o,-i,o - 5568//o,i,o,o,o - ^6m,o^3 - 1536//i,o,i,o,o) 
+ (X5 - 5/2x4 + 85/44x3 - 35/88x2 - 3/88xi)(-15048/5//i^2^) + - 5/2x4 + 1601/828x3 - 221/552x2 
-55/1656xi)(-3312//i,o,o6) + fe -5/2x4 + 8319/4292x3- 3497/8584x2 -265/8584xi)(-12876/57^o^2^) 
+ (X5 - 5/2x4 + 2559/1316x3 - 1097/2632x2 - 73/2632x0(1 1844^5) 

+ (X5 -5/2x4 + 71/36x3 - 11/24x2 - l/72xi)(20167^_i,o,-i,o,o) + (X5 -5/2x4 + 167/84x3 -27/56x2 

- l/168xi)(-1344i7o,i,o,i,o) + U5 -5/2x4 + 263/132x3 -43/88x2 - l/264xi)(52877_i,o,o^2) 
+ (X5 - 5/2x4 + 755/372x3 - 135/248x2 + 1 l/744xi)(-297677o,-i,-i,o,o) 

+ (X5 - 5/2x4 + 1063/516x3 - 203/344x2 + 31/1032xi)(412877o,-i,o,o,o) + - 5/2x4 + 25/12x3 
-5/8x2 + l/24xi)(-134477_i,o^3 + 115277_i,i,o^2 + 115277_i,i,o,o,o + 153677o,-i,o,i,o 

- 86477o,i,o6 + 1 152i7i,_i,o6 + 1 152i7i,_i,o,o,o + 1 152i7i,o,o,i,o - 1 152i7i,i,o^2 - 1 152i7i,i,o,o,o) 

+ (X5 - 5/2x4 + 379/180x3 - 79/120x2 + 19/360xi )(-2880//i,o,o,-i,o) + - 5/2x4 + 43/20x3 - 29/40x2 
+ 3/40xi)(-120//o,o,o,o,o) + fe - 5/2x4 + 545/252x3 - 125/168x2 + 41/504xi )(-2016//o,-i,o,-i,o) 
+ (X5 -5/2x4 + 13/6x3-3/4x2 + l/12xi)(-192//o,o^3) + (xj -5/2x4 + 485/204x3 - 145/136x2 
+ 77/408xi)(4087^o,o,o^2) -te -2522/1041x4 + 23633/12492x3 - 13051/24984x2 
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+ 443/8328xi)(24984/5^2^) + fe - 169/70x4 + 267/140x3-407/840x2 - l/120xi)(3360//o,i,o,o) 

+ (X5 -93/41x4 + 293/164x3 -805/984x2 + 295/984x1 )(3936//o,-i 6)- - 119/66x4+71/132x3 

+ 71/264x2 - l/264xi)(6336//i,o,o,o)-fe -98/59x4 + 63/236x3 +689/1416x2 - 131/1416xi)(5664//o,-i,o,o) 

-(X5 -477/322x4 -505/644x3 +2535/1288x2 -905/1288xi)(1288//o,o,o) + fe +49/10x4-769/60x3 

+ 751/120x2 +79/120xi)(480//_i,oAo) - (•^4 - 10x3 + 71/6x2 - 17/6xi)(96//i,o,o) + (x4 - 125/37x3 

+ 251/74x2 - 75/74xi)(2368//o,-i,o) + (^4 - 38/13x3 + 129/52x2 - 29/52xi)(-832i7o,i,o) 

- (X4 - 702/269x3 + 1297/538x2 - 431/538xi)(215277o,o) - (^4 - 90/37x3 + 139/74x2 

- 33/74xi)(355277_i^2) + (^4 - 367/152x3 + 1067/608x2 - 207/608xi)(4864^3) 

+ (X4 - 1 144/523x3 + 1249/1046x2 - 7/1046xi)(4184^2) + (^4 - 2x3 + 7/9x2 + 2/9xi)(-1728//_i^2ln2 

+ 57677-1,1,0,0 - 17287^1 ^2ln2 + 172877i,_i^2 + 576i/i,_i,o,o) - (x4 - 2x3 + 52/63x2 + 1 l/63xi)(2016i/_i,o,-i,o) 

+ (X4 - 2x3 + 7/8x2 + l/8xi)(1536//-i,o,i,o) + ix4 - 2x3 + 37/39x2 + 2/39xi)(4992//i,o,-i,o) 

+ (X4- 2x3 + 20/21x2 + l/21xi)(-1344//i,o,i,o) + (x4-2x3 + 129/131x2 +2/131xi)(6288//i ^3) 

+ (X4 - 2x3 + 94/93x2 - 1 /93xi )(-2976//_i,-i,o,o) + (X4 - 2x3 + 3 1 /30x2 - 1 /30xi )(- I92O//1 , 1 ,0,0) 

+ (X4- 2x3 + 25/24x2 -l/24xi)(10752//()_i_i,o) + (x4- 2x3 + 13/12x2 -l/12xi)(-2304//o -1,1,0 

- 384//o,i^2 - 2304//(),i,_i,o + 76SHqxiq) + (m - 2x3 + 74/63x2 - 1 l/63xi)(60487^_i,_i^2) 
+ (X4 - 2x3 + 5/3x2 - 2/3xi)(1728^2ln2 - 816//_i^3) 

+ (X4 + 20/7x3 - 271/42x2 + 109/42xi)(-6727/_i,o,o) + - 3/2x2 + l/2xi)(53767/_i _i,o - 2272i/_i,o 

- 1 152//_i,i,o + 607^0 - 192i?i^2 - 1 15277i,_i,o + 352i?i,o + 384i/i,i,o) + (x2 - xi)(20 - 321n^ 2 - 768Li4(l/2) 

- 384^2ln^2) + (-657/4^5 - 1 16^3 - 132^2 + 175/2^2^ - 205/2^2^ + 15677_i^2 - 327/1077_i^2^ + 12077_i,_i,o 

+ 88i7_i,o - 28i7_i,o^3 + 8677_i,o^2 + 4277_i,o,_i^2 - 14i/- 1,0,- 1,0,0 - 12/?-i,o,o - ^-1,0,0^2 

+ 1 1277_i,o,o,-i,o + 86i7_i,o,o,o - 1 1^-1,0,0,0,0 - 4877_ 1,0,0, 1,0 + 10077_ 1,0, 1,0,0 + 96i7_i,i,o + 24i7_i, 1,0^2 

+ 2477_i,i,o,o,o + 279/2i/o - 12977o6 + 309/2i7o6 + 159/477o^2^ + 6777o,-i6 + 3077o,-i^2 + 108i7o,-i^2ln2 

-234//o,-i,-i^2 -22//o,-i,-i,o,o - 12//o-i,o + 120//o_i_o^2 - 82//o_i_o_i_o + 6//o,_ 1,0,0 + 62//o_i_o_o_o 

+ 32//o,-i,o,i,o - 36//o,-i, 1,0,0 - 132//o,o - 8//o,o^3 - 75//o,o^2 - 123//o,o,-i^2 -432//o,o,-i,-i,o - 86//o,o,-i,o 

+ 259//o,o,- 1,0,0 + 144//o,o,-i,i,o + 245/2//o,o,o + 77/2//o,o,o^2 + 183//o,o,o,-i,o + 66//o,o,o,o - 9/2i/o,o,o,o,o 

- 74//o,o,o,i,o + 24//o,o,i^2 + 144//o,o,i,-i,o + 12//o,o,i,o -44//o,o, 1,0.0 - 48//o,o,i,i,o 

- 257//o,i^3 + 108//o,i^2ln2 - lOSHox-i^i - 36//o,i,_i,o,o - 64//o,i,o - 18//o,i,o^2 - 104//o,i,o,-i,o 

+ 98//o,i,o,o + 116//o,i,o,o,o + 4//o,i,o,i,o + 136//o,i,i,o,o + 16//i^2 + 513/10//i^2^ + 96//i,_i,o + 24//i,_i,o^2 
+ 24//i,_i,o,o,o - 176//i,o + 16//i,o^3 + 28//i,o6 - 186//i,o,-i^2 - 58//i,o,-i AO - 64//i,o,o + 55//i,o,o^2 

- 7677i,o,o,-i,o + 28//i,o,o,o + 71//i,o,o,o,o + 24//i,o,o,i,o + 32//i,o,i,o,o - 32//i,i,o - 24//i, 1,0^2 - 24//i,i, 0,0,0)^ W 
+ (-39^5 + 90^26 - 498/5^2^ - 138/5i/_i^2^ + 727/_i,o,-i^2 + 8/?_ 1,0,- 1,0,0 - 28//_i,o,o6 

+ 16//_i,o,o,-i,o - 12i?-i,o,o,o,o + 4Ho^3 - SHq^i + 307/5//o^2^ + 6877o,-i6 + 26477o,-i6 - 4877o,-i^2ln2 

- 168i7o,-i,-i^2 - 24i7o,- 1,-1, 0,0 + 647^0,-1,0^2 - 40/?o,-i, 0,-1,0 + 24i7o,-i,o,o + 4077o,- 1,0,0,0 + 16.ffo,-i,i,o,o 

- 12i7o,o^3 - 100^0,0^2 + 12i/o,o,-i^2 - 48i7o,o,-i,-i,o + 64i?o,o,-i,o + 3677o,o,-i,o,o - 877o,o,o + 6i?o,o,o^2 

+ 36/?o,o,o,-i,o -407^0,0,0,0 - 277o,o,o,o,o -4i7o,i^3 - 48//o,i^2ln2 + 48//o,i,-i^2 + 16i/o,i,-i,o,o - 16^o,i,o^2 
+ 16//o,i,o,-i,o + 8//o,i,o,o - 8/?o,i,o,o,o + 16//o,i, 1,0,0 + 138/5//i^2^ - 7277i,o,_i^2 - 8i?i,o,-i,o,o 
+ 28//i,o,o^2 - 16//i,o,o,-i,o + 12//i,o,o,o,o)<?(x)-' 

+ (-24^3 -8^2 -48^2ln2 + 269/5^2^ + 132//_i^3 + 384//_i^2 -48//_i^2ln2- 168//_i,_i^2 
-24//_i,_ 1,0,0 + 128//_i,o^2 - 168//_i,o,-i,o -48//_i,o,o + 168//_i,o,o,o + 128//_i,o,i,o + 16//-1, 1,0,0 
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- 60//o^3 - 224//o^2 + 164//o,-i^2 - 80//o,-i,-i,o + 288//o,-i,o + 108//o,-i,o,o + 96//o,-i,i,o + 48//(),o 
+ 122//o,o^2 + 228//o,o,-i,() - 180//o,o,() + lO/Zo.o.o.o - 144//o,o,i,o + 16//o,i^2 + 96//o,i,-i,o - 176//o,i,o 
-48//o,i,o,o - 32//o,i,i,o + 44//i^3 -48//i^2ln2 + 48//i,_i^2 + 16//i,-i,o,o + 112//i,o,-i,o 

+ 56//i,o,o - 56//i,o,(),o - 32//i,o,i,o + 48//i,i,o,oMx)^} 
+ en/Cpr^{(-16^3 - 296/9^2 + 16//-i^2 + 32//_i,-i,o - 592/9//-i,o - 16//-i,o,o + 1922/27//o + 4//o^2 

- 16//o,-i,o + 296/9//o,o + 8//o,o,o)^(^)} + <3(e^) ■ 

The one-loop part of Ti in Eq. (|20l) has first been given in Ref. O, and also T^j^* in Eq. (|22)) 
was known to (9(6) before yj. For the two-loop parts of Eqs. (|2T)) and (l23l) up to (9(6*^) we have 
found agreement with Ref. flU. All other terms, especially those at 0{e) in Eqs. (|2T|) and (|23l) . are 
new. 

In time-like kinematics above production threshold {Q^ > Am^) the conformal variable y is 
defined as 



, VQ^-VQ^^ (24) 



From Eqs. (|20l)-(|23l) the corresponding results for the form factor can be obtained by means of 
a suitable analytic continuation, x ^ y = -x, taking into account the (complex) continuation of 
Q^. This continuation is easily performed with the help of routines for HPLs [22] implemented in 
FORM [f38l . the only subtle point being logarithmic branch cuts, starting with 

\nx —> lny-l-i;r . (25) 

Thanks to Eq. (|25l) . both !Fi and Ti develop an imaginary part in the time-like region, 

r, = %ri+i3ri. (26) 

3.1 Asymptotic expansions 

Let us next study the asymptotic expansions of Eqs. (I20l)-(|23]). In the (space-like) high-energy 
limit, » m^, i.e. x — > 0, we abbreviate the logarithmically enhanced terms, 

L = l„g). (27) 

Then, employing algebraic properties of HPLs and keeping the scale /j.^ = we obtain from 
Eqs. (|20|)-(|231) the asymptotic expansions for the renormalized form factors. Here and in the 
following we set = 1/2. 

^i'^ = C^{^(2L-2)-L2 + 3L-4 + 2^2 + e(^Z.3-^L2 + (8-^2)Z^-8 + 2^2 + 46) (28) 

= C/{\{2L^-4L + 2) + -{-2L^ + SL^-{U-4^2)L + S-^2) + Il^-^L^ (29) 
I e 6 3 
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+ (f - 4^2)^2 - (y - 32^3)^. + 46 + 39^2 - 44^3 -48^2ln2 - -^^2' + ei^~L' + y 
/137 8 \ 3 /153 112 \ 2 /479 284 106 , 346 

- 24^2ln2 + 96^2ln^ 2- 160^2^ + 81n'^ 2 + 192Li4(l/2)- 12^2^3- 18^5)1 

f 1 / 11 11 \ 1//67 \ 49 \ 11 3 /233 \ . 

I 54 3 " ^ I 27 9^ 3 ^ " 5^ I 12 

/565 4 \ , /3337 11 \ 7 /70165 575 260 88 .\ 28745 

+ 6 - + —^2 - 26<-3 + + 6 6 6 

\ 54 3^ / I 54 2 ^ ^ / I 324 18 3 ^ 5 ^ / 
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- + + 12^2ln2 - 41n4 2 - 96Li4(l/2) - 48^2ln' 2 + ^^2' - 2^6 - 157^5) 

^ fl/2, 2\ 1/ 10, 10\ 2 , 19 . /209 4 \, 212 14 8 

+ Cr,n,{ — -L - - + L + — \--L^ + —L^ - + -6 + 6 - -6 

^ 'le2\3 3/ el 9 9 / 9 9 I 27 3^/ 27 9 ^ 3^ 



/1 . 4 47 , /281 \ . /5813 37 8 \, 2602 88 140 48 

+ e --L^ + + 6 ^2 - + — A - -6 L + <-2 (3 (2 ' 

l6 27 I 27 \ 162 9 " 3^ / 81 27^ 9 ^ 5 ^ 

f 2 . 19 , /265 \ 766 4 /I 4 47 , /281 \ , 

+ Cp — + — - + 2^2 + ^2 + 4 -L^ + { + ^2 

9 9 \ 27 ^ I 27 3^ \6 27 \ 27 ^ I 

/191 ^ 28 \ 2069 808 112 . ^ 92 4 ,\i , . 

^i'^ = ^C^{4L + e(-2L2 + 16L + 4^2) + e'(^Z.3-8L2 + (32-2^2)z.+ 16^2 + 8^3)} + 0(e^), (30) 

T^^ - ^C/|^(8L2-8Lj-8L^ + 34L2-(62-48^2)z. + 60^2-192^2ln2 + 16^3 (31) 
+ f(y - y + (74 - 16^2 - (249 - 1 16^2 - 192^3^1 

1 192 

- 20 + 492^2 - 1 152^2ln2 + 128^3 + 321n'* 2 + 768Li4(l/2) + 384^2 In^ 2 - ^-^2^ 

m2 f2 , 346 ^ 244 „^ . ^ „^ / 2 . /250 „ \ , 
+ 3^ + — ^ + 12 - —^2 + 96^2ln2 + 80^3 + ~\~^ ^^^ji' 

'8057 _ „ \, 1768 . 1504 

-384Li4(l/2) 

"i^^ f4,2 100, 8 / 4 , 148 ,2 /1922 ^ \, 296 

+ — rCp?i, -l? L - -6 + e — + 1} - + 44-2 L ^2 - 16^3 

Q2 F g 3" I 3 9 I 27 ^ I 9 ^ ^ 

n? (A 2 100 136 „ / 4 , 148 , /1922 ^ \ 764 128 
g2 f\3 9 3 ^ \ 3 9 \ 27 ^ I 9 9 ^ 
128 



+ (-^ + 38^2 - 48^3 )l + 78 - -^^^2 + 576^2 In 2 - 264^3 - 192^2ln^ 2 + -^^2^ - 161n'* 2 



+ 64^2ln2-— ^3)} + C(e') 



For the magnetic form factor only the first power in irP-/Q^ has been kept in the expansion. 
In Eqs. (|29|) and (|3TI) . we have also separated the contribution proportional to C^T^ coming from 
the heavy-quark loop in diagram QL. This part is finite after renormalization and proportional to 
a single power of Cp only, thus it enters at 0{e^) and O(e^) at two loops. Moreover, due to the 
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diagram 
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8 


result for T\ 


result for T'l 


L 




{Cr-'-tf 


Eq. dill) 


Eq.(l35l) 


C 




[Cf-^-YCf-Ca) 


Eq. 


Eq.dlTl) 


Qv 


(Cf-^)Cf 

\ ^ f 




Eq. dSH]) 


Eq.dH 


s 


F 




Eq. dlOl) 


Eq.dlB 


QL 


Cf 


Cf-^ 


Eq. 


Eq.dll) 


GL 


CfCa 


[Cf-^)Ca 


Eq. dUl) 


Eq. d45l) 


GV 


\CpCa 




Eq. (1461) 


Eq.dlH) 



Table 2: The color coefficients multiplying the result of the individual Feynman diagrams of the quark form 
factor (cf. Fig. ^ for a current in the singlet 1 or octet 8 representation of the color-SU(A'^c) gauge group. 
The scalar functions are given in the respective equations. 



magnetic form factor ^Fj^^ in Eq. (l30l) being finite also the terms proportional to the number of 
light quarks are finite in Eq. (|3TI) . 

The expansions d28l) and (l29l) for Ti in combination with the exponentiation of the heavy-quark 
form factor provide an independent consistency check on the correctness of our result, in particular 
of the new terms of order e at two loops in Eq. d2T|) . A detailed discussion of this aspect will be 
presented in the next Section. 

Let us instead now turn to the limit ~ 4m^, i.e. j — > 1, which provides the threshold expan- 
sion in time-like kinematics, cf. Eq. d24l) . The relevant small parameter in which we expand here 
is the heavy-quark velocity jS, 




^=a|1-^- (32) 

In this particular limit, i.e. ;0 ~ 0, the well-known Coulomb singularities appear. The inverse powers 
of contributing to the form factor have already been directly addressed in investigations of non- 
relativistic QCD (e.g. [iBQll^Ol ). Without repeating too much of the discussion in the literature, it is 
perhaps instructive, to study the anatomy of Coulomb two-loop corrections based on the individual 
diagrams which we have at our disposal. To that end, we list below the individual results for the 
threshold expansion of the bare Feynman diagrams from Fig. [2lup to terms of Oifi) and terms of 
0{e) as they enter in the computation of the bare form factor Tf'^^\ 

<^b,(2) ^ 2tI^^ + Tl^^KTf^^ + 2TI^^^ + 2Tf^^ +Tf^ + tI^^ (33) 

with / = 1,2 and the symmetry factor of 2 has been added for the diagrams GV, QV and S . It is 
also understood, that the respective color factor multiplying each diagram is taken from column 1 
in Tab. [21 i.e. for a color-singlet current: 

<jf^ = 1{JJ_^,,9_3 3 1| 1|1 /J X M 
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Due to the time-like kinematics we naturally encounter imaginary parts. We should point out 
again, that the diagram QL contains both light quark loops and one heavy quark loop, while GL 
is the gluon and the ghost loop contribution to the gluon self-energy. 

While our calculation (|33] ) has been performed for a color-singlet (1) current coupling within 
the color-SU(A^t) gauge group, nevertheless the results of Eqs. (I34l)-(|471) can also be applied to 
considerations of color-octet (8) currents. In the latter case the relevant color factors are listed in 
column 8 of Tab. [21 Such decompositions are needed, for instance, in the description of heavy 
quarkonium states (see e.g. [411 |42l |43]|). At one-loop, i.e. for the diagrams OneV and OneS in 
Fig. [T] the transition from a quark-pair in a color-singlet state to one in a color-octet state amounts 
to the simple replacement of Cf {Cf - Ca/2). At two loops, this is not true for the individual 
Feynman diagrams, see Tab. [2l However, for the complete two-loop Coulomb corrections again 
the same simple replacement rule holds, as we want to illustrate next. 

To that end, we present our results in the form as they enter in physical processes, 

X ^ (qq)i, 1 = 1,8, (48) 

where {qq)i represents a (heavy) color-singlet or octet final state. Near threshold, considering only 
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Coulomb enhanced contributions, the relevant Bom cross section cTj*^^ factorizes, so that 

crx^igcik = <Tf^(l+a,A['KalAP), (49) 
where the Aj^ can be expressed in terms of the threshold expansions of the renormalized form 



factors above threshold, cf. Eq. (1261) . For the color-singlet case (,qq)i we have from Eq. (|33l) . 

aJ^^ = 2(%ri^'^ + %r^^^), (50) 
Af = (%ri'^+%r^'^f+(3rl'^ + dr^^f+2(%ri^^ + %r^^). (5i) 



Upon adding Eqs. (I34l)-(|47]) with the appropriate color coefficients from Tab. [2] and performing the 
necessary renormalization, i.e. the threshold expansion of the time-like counter-terms (fT2l) -(fT3l). 
we find the Coulomb corrections for the color-singlet final state as 

(52) 

31 10 2 n"^] 

-Q -16C,- -n, - 2^oln 2 - 2ySo In^j - + -C,-^ . (53) 

Note, that all divergencies in e as well as all higher inverse powers in yS have canceled in Eqs. (I52l) 
and (l53l) . so that we retain at most at L-loops. For the octet final state (qq)8 we can 

evaluate Eq. (l33l) with the color factors appropriately replaced according to Tab. [21 The resulting 
expressions for A^ * and Ag^^ after these modifications are then indeed obtained as in Eqs. (l52l) . 
(l53l) with the replacement of the color factor 2C^ (2C^ ~^a^- ^^^^ complete agreement 
with results from potential non-relativistic QCD (e.g. |[39ll40l ). 



Ai^) = 


= 2C/- 


Af) = 


- 2C^|(: 



4 Applications 

In this Section we present several applications of the new results for the heavy-quark form factor. 
Particular emphasis will be put on new three-loop predictions, because the order e results at two 
loops are necessary ingredients for a full analytic three-loop calculation. The applications will 
mostly be concerned with the electric form factor !Fi . We discuss the property of exponentiation 
as well as its role in the determination of massive n-parton gauge amplitudes in the high-energy 
limit from massless ones. 



4.1 Exponentiation of form factor 

Let us start by recalling briefly the exponentiation for the heavy-quark form factor. This feature is 
based on the universality of soft and coUinear radiation and the respective singular terms in T'l, i.e. 
the poles in e and the large logarithms ln(m) of Sudakov type [4], so that fi fulfills the following 
evolution equation (see e.g. [i5ii6il), 

2 d ^ ^IQ' \ 1 Jm^ \ l/Q^ \ 
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lim// —K\—,as,e\ = n —G\—,as,e\ =A(a^). (55) 



where the dependence on the various scales has been made exphcit. QCD factorization at the scale 
H allows to separate in Eq. (|54|) the dependence on the hard scale Q from that on the heavy-quark 
mass m. To logarithmic accuracy, the former rests entirely in the function G while the latter is 
associated with the function K. Both functions, G and K, are subject to renormalization group 
equations [[51 im governed by the same (well-known) cusp anomalous dimension A 11451 146ll47ll . 

ri TiKy —,as,e\ = n — — 

m-»0 af/i^ \//^ / rf/i^ \yU^ 

In dimensional regularization, the solution of the evolution equation (l54l) proceeds in complete 
analogy to the construction for the form factors of massless quarks and gluons, see for example IfTSl 
|48l|49l|50l. In the massive case, the necessary integration constant is parameterized through the 
matching function C. For the (ultraviolet) renormalized massive form factor T\ with space-like 
virtuality = -Q^ < and in terms of the renormalized coupling Usip^) we arrive at 

ri{^,'^,asiM\^ = Cia(M\e),e) (56) 
i / jiGmfi\e)) + K(a(^^W/Q\e))- J ^A(a(V,6)) 



xexp 



where all quantities on the right hand side are defined in D dimensions (note a change of sign in 
the definition of K compared to Ref. [6J). They are functions of the D-dimensional strong coupling 
a and can be expressed in terms of the ordinary coupling as in four dimensions (see e.g. [il5il ) like 
the functions C(a), G{a) and K{a). 

The coefficients in the perturbative expansion of the cusp anomalous dimension A and the 
function G agree with those of the form factor for massless quarks, see e.g. [fT5l |49ll . Due to 
the heavy-quark mass acting as an additional regulator in the coUinear limit the coefficients of the 
infrared counter term K and also the matching function C take, on the other hand, particular values. 
Their respective coefficients read 

K^^^ = -2Cp, (57) 
= c/(-3 + 24^2 - 48^3) -C^C^(^ + 30^2- 60^3) + «/Cf(^+ 4^2), (58) 

= Cp{4 + ^2) + eCp{s + ^^2-^^3) + e'Cpl\6 + 2^2-l^3 + ^^2') + 0{e'), (59) 

n^ 7(241 163 , / 557 399 286 

= C/{— + 26^2-6^3-48^2ln2-— ^2' + e(— ^ + — ^2 + — ^3-24^2ln2 (60) 

A J 1363 , 146 

+ 192Li4(l/2) + m^2 + 96^2\r^2- —^2^ - —^2^3 - 6^5 

f 12877 323 89 . ^ 47 , /629821 7805 1162 209 , 

83 

- 96Li4(l/2) - 41n'* 2 + 12^2ln2 - 48^2ln^ 2 + y ^2^3 - 208^5 j j 

f 1541 37 26 / 46205 673 182 49 2\ 

+ n,CA ^2 6 + q ^7 ^3 CiW + Oie^] . 

' 324 9 ^ 9 ^ I 1944 54 27 ^ 5 ^ // ^ ' 

Eqs. (|57l)-(|59l) have been obtained before The 0{e) term of C*^^ in Eq. (|60l) is a new result. 
Note also that we have changed the definition of C*^^^ compared to Ref. [^, because matching at 
two loops has now been extended to 0{e) consistently. 
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With these values for A, G, K and C and, upon performing the integrations over ^ and A in 
Eq. (|56l) the exponential for !Fi generates all large logarithms in the heavy-quark mass m (even 
including higher orders of e). In particular, we are able to reproduce the high-energy expansions 
through two loops in Eqs. (|28l) and (|29l) . This result provides an independent check of our result, 
because the exponentiation (l56l) relies on functions A and G that have been derived in massless 
computations. 

Let us at this point make two remarks on the exponential (l56l) . First, the derivation of Eq. (I56l) 
relies on the heavy-quark masses (acting as infrared regulator) to be attached to external lines. By 
construction, we have omitted all contributions from diagram QL in Fig. [21 with a virtual heavy- 
quark in the loop in Eq. (|56l) and the comparison to Eqs. (|28l) . (|29l ). These contributions are finite 
after performing the ultraviolet renormalization, but they still do contain Sudakov logarithms L, 
see Eq. (ITT]) . It is well-known, that these remaining logarithms also obey an exponentiation similar 
to Eq. (|56l ). see e.g. Ref. (61 and related discussions in electro-weak theories |l5Tl[52l|. Here, we do 
not elaborate on this point further. 

Second, in the exponential expression Eq. (|56l ) for the form factor we have used the standard 
MS coupling running with light flavors. In order to compare Eq. (|56l) or rather its expanded ver- 
sion to the fixed-order results Eqs. (I20l) and (|2T]) which also employ the MS-scheme, but a running 
coupling with a total number of flavors = n^ + l, one has to apply the decoupling relations [l53ll . 

The necessary decoupling constant for a? at flavor thresholds is known to 0(al) ||54l|55]| (see also 
Refs. [[561I571 for C)(Q'j)). To relate the two results, we use the following relation for as, 

„r/'.(„r''f{4L,-.(iL^ift)-e^(iL^if.i,-2^3)) (61) 

/ (n,.)\3^4 10 16 15 

+ [c^s ' ) { - -CaL, - 2CpL, + -C^ - -C, 

/4 . 10 9 32 29 4 86 5 35 \ 



a/ 



■ e 



where is the standard MS coupling for quark flavors expanded in terms of for = -I- 1 
flavors, both evaluated at the scale jjp-, and 

Z,, = l„(g. (62) 

Eq. (|6TI) uses the pole-mass m and the higher order terms in e in Eq. (|6TI) are needed for the higher 
order expansion coefficients in Eqs. (|28|) and (|29|) . 

Now we are ready to employ Eq. (l56l) to obtain a partial new three-loop prediction for T\ in 
the limit :^n^ . With the necessary expansion coefficients for A and G to third order, we obtain 
at the scale fj. = m, 

= C/{^(^l3-4L2 + 4L-1) + 1(-2L4 + 10L3-(22-4^2)Z.' + (22-8^2)Z.-8 + 4^2) (63) 
^ ~ T^' ^ ("T ~ ^^'K ~ (^^ " ^^^'K ^ (^^^ ^ ~ -96^2ln2 - ^^2')z. 
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fl/ 4 4 \ (20 8 \ 1 

6e 

As detailed above, we have omitted the finite contributions arising from virtual heavy-flavor lines, 
i.e. diagram QL in Fig.[2]with the heavy-quark in the loop and higher order generalizations thereof. 
All singularities in e however are controlled by the decoupling relation (|6TI) which also gives rise 
to the terms proportional to C^^,C^C^ and C^n^. The latter are explicitly displayed in Eq. ([63]) . 

The accuracy of the prediction for T^^^ in Eq. (|63] ) is limited by lacking knowledge about 
the three-loop function K^^^ controlling the infrared sector of the massive theory. K^^'' has been 
kept explicitly in Eq. (|63l) above. The constant terms of order 0{e'^LP) in T'^^ contain the sum of 
G*-^-* including its order (9(6)-part which became available recently Il58l . as well as the currently 
unknown three-loop matching function C^^^ . Therefore we have truncated T'^'' as given in Eq. (|63l) . 

4.2 Massive w-parton amplitudes 

Let us again briefly summarize the key features of relating loop amplitudes with massive partons 
to massless ones. QCD factorization implies that a massive amplitude A1^'"^ for any given phys- 
ical process shares essential properties with the corresponding massless amplitude Al*^™"^^ in the 
limit when all kinematical invariants are large compared to the heavy-quark mass m. One typi- 
cally encounters in M^"^'^^ two types of singularities, soft and coUinear ones, which are related to 
the emission of gluons with vanishing energy and to coUinear parton radiation off massless hard 
partons, respectively. In dimensional regularization, these appear explicitly as factorizing poles 
in 6, see BH HOl, and the associated anomalous dimensions are currently subject of active stud- 
ies [|50l |59l [601 [6T1 [621 [63ll . In the massive case, the soft singularities remain in A1*^™^ as single 
poles in e and the heavy-quark mass m screens some of the coUinear singularities giving rise to 
logarithms in m, see [^,T[. Recent progress on the mass dependence in the relevant soft anomalous 
dimensions has been achieved in Refs. [il3l[T4l[64l . 

The common physical origin of the singularities lends itself to a proportionality between At^'") 
and M^"^^^^ so that QCD factorization provides the remarkably simple and suggestive direct rela- 
tion M 

M^'") = W (zf;;'l«^)^xM('"=o). (64) 

ie (all legs} 

The universal multiplicative factor Z^"^^^^ is process independent and depends only on the external 
parton. Then, Z^*"'*^^ is defined simply as the ratio of the on-shell heavy-parton form factor and the 
corresponding massless on-shell one. E.g. for an external heavy quark q, we find 

where only the electric form factor !Fi enters in the high-energy limit 2^ » m^, since Ti vanishes 
for massless quarks. Here Us is evaluated at the scale fj^. The process independence of Z^'"''^^ is 
manifest, because the hard scale drops out in Eq. ([651) leaving only a function of the ratio of 
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scales fj.^/m^. The definition (1651) of Z^'"'^\ however, excludes terms with explicit dependence on 
the number of heavy quarks (e.g. the virtual heavy loop in diagram QL in Fig. |2l). At two-loops 
(and beyond) one needs additional process dependent terms for their description [[8]|. 

Applications of the formalism so far include the computation of the (virtual) two-loop QCD 
corrections to hadronic top-quark pair-production [fTTl [T2l in the limit when all Mandelstam in- 
variants are large, 5, » fn^- In the same approximation, the derivation of the two-loop 
QED corrections to Bhabha scattering for a small electron mass lUg has been performed [8J (see 
also ll65l [66l [671 for earlier work). These applications have also covered the complete depen- 
dence on virtual contributions for heavy fermions in loops, mostly by means of direct calcula- 
tion UMl EH EOl Ell IH . 

With the help of Eqs. (1631) . (I65l) and the known results (15] of the poles for the massless quark 
form factor we can extend the perturbative predictions for Z*^*"'*^^ to higher orders in a 5. Defining 

,/=i j=i 

we can derive the three-loop coefficients Z^^^ up to the single poles in e. 

Keeping the currently unknown function K^^\ see the exponential in Eq. (I56l) . we obtain 
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+ 161n42 + 384Li4(l/2) - ^6' - 8666 - y + C/C^I^^L^ + (-^ - 96)^;^ 
1157 899 \ 2 /4507 521 , , 140 297 2\ 138403 



/1 15/ 899 \ 2 /450/ 521 . „ 140 29/ -,\ 

+ + <2 - 1356 U^r, + + 6 + 14461n2 + 6 6 + ■ 

I 12 12 ^ "/ \ 27 4 ^ ^ 9^5^/'^ 486 

19211 .„ . _ 2441 „, . . 1633 



+ ^^6 + 4861n2 - —6 - 9661n' 2 + — 6^ - 8111^2 - 192Li4(l/2) + 4166-4166) 
2/245 536 44 88 2\ 2/ 107 , 17 . /73 67 \ . 

I 27 2 ^ 9 ^ / ''^ 243 54 ^ 9 ^ 20 ^ / 
/418 80 56 \ 28 1 ni1 . n 

where is defined in Eq. (|62l) . Eq. (l67l) has been obtained from an all-order in € formula in 
Ref. [|6|| based on Z^^^j being identical to the virtual contribution (prior to coUinear factorization) to 
the perturbative fragmentation function of a heavy quark [[73J . 

According to the discussion above all contributions arising from virtual heavy-flavor loops 
are omitted here and it requires substantially more effort to determine them in massive scattering 
amplitudes, see e.g. Refs. [fm[T2l for the two-loop QCD corrections to hadronic top-quark pair- 
production. It should be clear though, that the highest power in = + is always governed by 
the renormalization of the coupling constant a?. Thus, terms proportional to the maximal power 

can be predicted on the basis of the Z'^'"''^^ -factor for massive scattering amplitudes as well. 
However, sub-leading terms in remain inaccessible or require extending the results of Ref. dH. 

With the recent progress in the computation of the massless quark form factor [158117411751 . one 
can hope to extend Eq. (|69l) to the constant terms as well. However, to that end, we would need to 



33 



know also the high energy-limit of T'^ even to order 0{e^) as well as that of to order 0{e^). 
In terms of the functions K and C governing the exponentiation (l56l) this requires besides K^^'' 
also the two- and three-loop matching coefficients C^^^ and C^^^ to sufficient depth in e. Of these 
quantities, at least the coefficient C'^^-' is not entirely out of reach, because the necessary two-loop 
master integrals are known to sufficient depth in e (see Sec. 111). 

In closing, we comment once more on the importance of Eqs. (l64l) and (I67l)-(l69l). We can 
derive all logarithmically enhanced terms and even the constant, i.e. mass -independent terms of 
A1*^'"\ provided Al*^'""^-* and all Z^^^^ are known to sufficient depth in e. For massless scattering 
amplitudes we do understand the general structure of singularities [|9l[l0||. This allows to predict the 
pole structure in e of massless amplitudes at any order based on a small number of perturbatively 
calculable anomalous dimensions (see also [|50l [611 [62l [631). This is a constructive approach to 
y^(m=0) Thus, with the help of Eq. dill), the singular limit of M^"'=^^ and the factors Z^'^ ,i=l,...,3 
from Eqs. (I67l)-(|69l) one can study the logarithmically enhanced terms in an amplitude Al*^*"^ with 
external massive fermions. It is now, for instance, a straightforward exercise to investigate explicit 
results for the virtual contributions to massive-fermion scattering amplitudes in the small-mass 
limit at three loops, an obvious example being e.g. the three-loop virtual QED amplitude of Bhabha 
scattering to logarithmic accuracy. We leave these issues to future work. 

5 Conclusions 

The form factor of heavy quarks has been studied. Its radiative corrections in QCD are of great 
relevance for precision predictions of many observables at colliders. Moreover, on the formal side, 
they exhibit very interesting structures of quantum field theory. 

We have presented a fully rigorous calculation of the two-loop QCD corrections to the electric 
and magnetic form factors !Fi and 9^2 in dimensional regularization. This allows for an independent 
check of the previous computation reported in the literature In addition we have obtained new 
two-loop results at order e. To that end, we have determined the Laurent series of the necessary 
master integrals to very high order in e. The expansion of our result in various limits has delivered 
new insight into the anatomy of the two-loop form factor near threshold and the Coulomb singu- 
larities for color-singlet and octet currents. The high-energy limit on the other hand provided new 
tests of the exponentiation in D dimensions as well as an improved three-loop prediction for the 
form factor to logarithmic accuracy in the heavy-quark mass m. We have also discussed the conse- 
quences for massive scattering amplitudes, where we have been able to advance towards three-loop 
predictions for massive n-parton amplitudes in the small-mass limit based on the known singular- 
ity structure of the massless result. As the heavy-quark form factor and its QCD corrections are 
essential ingredients in high precision theory predictions, we believe, there will be several other 
future applications, where our results will prove to be useful. 

Files of our results can be obtained from the preprint server http: //arXiv.org by down- 
loading the source. Furthermore they are available at f3T| or from the authors upon request. 
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